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EVALUATION 


The  objective  of  this  effort  is  to  develop  and  validate  a  simulation 
of  the  space-based  radar  (SBR)  phased  array  lens  system.  This  report 
contains  the  simulation  theory  and  implarentation  methods  used  to  set  up 
the  simulation.  The  key  accomplishment  of  this  effort  has  been  to 
establish  a  new  method  to  perform  an  analysis,  to  the  accuracy  desired, 
of  a  large  SBR  phased  array  lens  system  that  includes  all  electrical 
factors  that  influence  the  RF  far  field  pattern.  Also,  the  key  parts  of 
the  simulation  have  been  validated. 

This  effort  is  part  of  RADC  TP0-3A. 

(R.ajrf/f. 

RUSSELL  C.  STEENROD,  Capt  USAF 
Project  Engineer 


SECTION  1 


INTRODUCTION 


The  objective  of  this  contract  is  to  develop  and  validate  a  computer  program 
for  simulation  of  space-based  radar  (SBR)  phased  array  lenses.  This  simulator  will 
be  used  to  analyze  the  performance  of  candidate  lenses.  The  parameters  of  primary 
concern  are  radiation  power  gain  and  pattern  characteristics. 

A  typical  SBR  is  shown  in  Figure  1-1.  The  transmitter  provides  the  space  feed 
that  illuminates  the  lens  during  radar  transmit.  The  space  feed  may  be  in  the  form 
of  several  independently  controlled  beams  that  provide  lens  radiation  pattern  shap¬ 
ing,  adaptive  nulling,  and  time  delayed  lens  sector  illumination.  (The  time  delayed 
compensation  is  required  because  of  significantly  different  transmitter-to-lens  path 
lengths  to  different  points  on  the  lens).  Beam  steering  and,  possibly,  amplifica¬ 
tion  are  performed  by  the  lens. 

A  lens  comprising  two  arrays  of  metallic  radiators  sandwiching  a  ground  screen 
currently  is  being  modeled  in  the  simulator.  A  typical  lens  of  this  type  is  shown 
in  Figure  1-2;  for  simplicity,  dipole  radiating  elements  are  indicated.  The  simula¬ 
tor  also  will  be  capable  of  analyzing  other  element  types;  e.g.,  folded  dipole, 
bowtie,  turnstile,  and  parasitic  (Figure  1-3).  Focusing  and  scanning  of  the  lens- 
transmitted  main  beam  is  accomplished  by  electronic  modules  interconnecting  the 
radiating  elements  between  the  illuminated  and  nonilluminated  arrays  (Figure  1-4), 

These  modules  also  may  provide  power  amplification. 

An  accurate  analysis  of  a  large  space-based  phased  array  lens,  comprising  hun¬ 
dreds  of  thousands  of  elements,  is  a  formidable  problem.  Infinite  periodic  array- 
theory  alone  is  not  adequate  when,  for  example,  edge  effects  must  be  considered. 

A  space-based  lens  is  expected  to  have  stringent  constructural  and  deployment  con¬ 
straints  that  will  result  in,  at  best,  a  piecewise  periodic  lattice;  e.g.,  a  lens 
composed  of  wedge-shaped  gores  may  have  nonperiodic  spacings  in  the  vicinity  of 
gore  interfaces.  Another  source  of  nonperiodicity  may  be  inherent  in  the  electronic 
modules.  They  may  not  be  identical  in  the  event  of  selected  power  failures. 

(Simple  progressive  phasing  between  module  settings,  however,  can  be  treated  with 
a  variation  of  ''standard''  infinite  array  theory  as  described  in  this  report.) 

The  large  number  of  elements  in  the  lens  arrays  and  the  expected  nominal  peri¬ 
odicity  within  sections  of  each  array  suggests  obtaining  a  first  approximation  of 
element  currents  from  an  infinite  array  analysis  of  each  section  since,  within  a  ; 
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section,  periodic  elements  and  nearly  uniform  module  amplitude  settings  may  be 
assumed.  It  is  expected  that,  for  limited  scan  angles,  the  currents  on  many  ele¬ 
ments  within  a  section  will  assume  these  first-order  solutions.  The  phased  array 
simulator,  therefore,  makes  extensive  use  of  infinite  array  theory. 

A  successive  approximation  method  is  employed  whereby  infinite  array  analysis 
provides  a  first  approximation  of  array  currents  and  patterns.  Higher  approxima¬ 
tions  then  are  obtained,  with  relative  ease,  by  methods  that  draw  upon  these  first 
approximation  currents.  The  higher  approximations  provide  corrections  to  the  cur¬ 
rents  on  elements  in  the  vicinity  of  discontinuities  in  module  amplitude  settings 
or  periodicities;  e.g.,  lens  edges  and  section  interfaces.  Scattering  from  nearby 
structures,  such  as  supports  and  hinges,  also  may  be  accounted  for  in  the  higher 
approximations. 

The  first  approximation  patterns,  although  involving  large  numbers  of  elements, 
are  computed  from  closed  form  expressions.  Since  the  number  of  elements  with  cur¬ 
rents  perturbed  from  their  infinite  array  values  is  expected  to  be  relatively  small, 
higher  approximation  "correction  patterns"  may  be  computed  by  direct  summation  of 
the  field  from  each  element.  The  superposition  of  the  first  approximation  and  cor¬ 
rection  patterns  results  in  the  overall  pattern.  This  approach  takes  maximum  advan¬ 
tage  of  the  nominal  periodicities  and  largeness  of  the  arrays.  It  is  as  computa¬ 
tionally  manageable  in  the  radar  receive  mode  as  in  the  radar  transmit  node  without 
resorting  to  reciprocity.  This  is  important  since  the  modules  may  well  be  nonrecip¬ 
rocal.  This  approach  also  may  avoid  the  need  for  a  fast  Fourier  transform  (FFT) 
for  pattern  computation.  An  FFT  is  not  practical  if,  for  example,  fine  angular 
pattern  resolution  in  the  vicinity  of  a  null  is  required. 

An  infinite  array  analysis  usually  assumes  plane  wave  excitation.  During  radar 
transmit,  the  illuminating  field  is  expected  to  differ  smoothly  from  a  plane  wave 
across  the  face  of  the  lens.  This  difference  is  taken  into  account  by  expanding 
the  beamformer  field  in  plane  waves  and  analyzing  the  lens  separately  for  each 
plane  wave  component. 

The  first  approximation  infinite  array  method  employed  here  is  an  extension  of 
a  moment  method,  plane  wave  expansion  technique  recently  developed  by  Munk  and 
Burrell.1  It  is  applicable  to  lens  arrays  composed  of  arbitrarily-shaped  wire  radi¬ 
ating  elements.  The  elements  may  be  inclined  with  respect  to  the  array  plane.  Feed 
line  scattering,  amplitude  and  phase  adjust  modules  between  the  arrays,  and  an 
imperfect  ground  screen  also  are  accounted  for.  In  particular,  this  technique 


permits  the  modules  to  be  progressively  phased,  although  the  module-to-module  ampli¬ 
tude  adjustment  is  assumed  uniform  throughout.  The  plane  wave  expansion  technique 
facilitates  computing  the  array-to-array  coupling  through  the  imperfect  ground  screen 
and  suggests  straightforward  extensions,  not  discussed  here,  applicable  to  dielec¬ 
tric  support  sheets  (e.g.,  Kevlar  and  Kapton)  on  which  the  radiators  may  be  mounted. 

The  infinite  array  lens  analysis  facilitates  computation  by  providing  the  first 
approximation  solution  to  a  perturbation  problem;  it  also  provides  clues  to  the  per¬ 
formance  of  several  important  finite  (but  large)  lens  characteristics.  For  example, 
the  array  active  impedance  variation  with  module  phase  setting,  an  important  para¬ 
meter  for  module  designers,  can  be  assessed  directly.  The  effect  of  imperfect 
ground  screen  "feedback"  from  the  target  side  to  the  feed  side  of  the  lens  also  can 
be  observed. 

A  general  discussion  of  the  simulator  and  its  current  design  status  is  contained 
in  Section  2.  A  detailed  theoretical  discussion  follows  in  Section  3.  The  program¬ 
ming  and  data  flow  are  discussed  in  Section  4.  Preliminary  results  and  validation 
are  given  in  Section  5. 
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Figure  1-4.  Cross  Section  of  Lens 


SECTION  2 


DISCUSSION 


2.1  SIMULATOR  CAPABILITIES 

The  simulator  for  analyzing  candidate  SBR  phased  array  lenses  must  predict  the 
RF  impact  of  the  lens  on  the  following  SBR  gain  performance  factors: 

1.  Beam  scan  coverage, 

2.  Main  beam  and  sidelobe  beam  widths,  positions,  and  levels, 

3.  Null  depths  and  positions, 

4.  Effective  radiated  power, 

5.  Available  received  power,  and 

6.  Bandwidth 

These  performance  factors  must  be  evaluated  as  a  function  of  several  lens  parameters: 

1.  Structure  fabrication  and  deployment  variations 

a.  Array  lattice  tolerance 

b.  Element  orientations  and  locations  relative  to  ground  screen 

2.  Space  environment  effects  (heating,  etc.);  e.g.,  lens  warping 

3.  Module,  phase  and  amplitude  variations 

a.  Production  tolerances 

b.  Failures 

The  sumulator,  besides  analyzing  large  periodic  planar  arrays,  thus  must  account  for 
the  effects  of  lens  warping,  discontinuous  periodicities,  and  edge  elements. 

The  ground  screen  for  a  candidate  lens  is  expected  to  be  of  coarse  mesh  construe 
tion  to  minimize  overall  weight  and  temperature  gradient  effects.  For  an  active  lens 
where  the  nodules  impart  gain,  the  isolation  between  illuminated  and  nonilluminated 
arrays  may  not  sufficiently  suppress  the  ’’feedback"  penetration  coupling  between 
these  arrays.  This  interarray  radiation  coupling  must  be  included  in  simulations. 

The  array  lattices  may  be  triangular  or  rectangular;  thus,  simulations  must 
apply  to  general  rectilinear  arrays.  The  array  elements  may  be  turnstiles,  dipoles, 
folded  dipoles,  etc.  They  may  have  parts  inclined  with  respect  to  the  array  plane; 
e.g.,  V-type  dipoles.  The  effects  of  scattering  from  the  feed  lines  and  unbalanced 
mode  feed  line  radiation  also  must  be  considered. 

A  phased  array  lens  simulator,  therefore,  must  consider  many  effects  in  accu¬ 
rately  modeling  SBR  lenses.  This  capability  of  the  simulator  leads  naturally  to  its 
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OPTIMIZED  RADIATOR  PHYSICAL  CONSTRAINTS 

RADIATOR  TYPE  GROUND  SCREEN  CHARACTERISTICS 


Figure  2-1.  Simulator  as  a  Elusion  Tool 


use  as  a  lens  design  tool  as  uell  as  a  means  for  evaluating  designs.  The  interac¬ 
tion  between  the  simulator  and  the  designing  membranes  (radiating  elements  and 
ground  screen),  modules,  and  space  feeds  is  shown  in  Figure  2-1.  The  membrane 
designers  can  decide  on  array  lattice  size  and  type  (rectangular,  equilateral 
triangular,  etc),  ground  screen  transmission  characteristics,  and  radiating  ele¬ 
ment  physical  constraints.  The  simulator  then  can  be  used  to  "optimize"  the  choice 
of  a  radiating  element  type,  shape,  and  size  within  these  constraints  such  that 
active  impedance  variation  with  scan  is  minimized.  The  latter  is  important  in 
module  design. 

2.2  CURRENT  STATUS,  FUTURE  IMPROVEMENTS,  AND  EXTENSIONS 

The  status  of  key  components  of  the  phased  array  lens  simulator  is  shown  in 
Table  2-1. 

The  theory  category  is  subdivided  into  "general”  and  "detailed.”  A  detailed 
description  includes  all  equations  down  to  the  lowest  level.  A  completed  "code" 
refers  to  flow  charts  as  well  as  a  computer  program.  A  satisfactory  "debug”  indi¬ 
cates  the  code  is  implemented  and  providing  reasonable  results.  A  successful 
"validation"  connotes  agreement  with  other  methods  in  solving  special  problems. 

These  other  methods  could  be  theoretical,  computational,  or  experimental.  Validated 
results  are  presented  in  Section  5.  A  detailed  plan  for  validating  all  aspects  of 
the  simulator  is  being  prepared  and  will  be  published  in  a  separate  document. 

The  only  component  of  the  simulator  that  is  not  expected  to  be  completed  under 
the  current  contract  is  "support  structure  scattering."  A  major  effrot  is  required 
here  to  implement  a  geometrical  theory  of  diffraction  within  the  higher  approxima¬ 
tion  moment  method  code.  Another  extension  to  the  simulator,  not  currently  envi¬ 
sioned,  that  fits  well  within  the  general  framework  of  a  plane  wave  expansion  moment 
method,  is  the  capability  of  including  the  effects  of  dielectric  support  sheets  such 
as  Kevlar  or  Kapton.  A  number  of  improvements  for  increasing  the  computational  effi¬ 
ciency  and  data  handling  also  can  be  investigated  in  subsequent  efforts;  e.g., 
optimizing  the  choice  of  current  expansion  functions  in  the  infinite  array  moment 
method  technique.  Other  methods  for  improving  efficiency  will  be  implemented  if 
time  permits  (Section  3.2."). 
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SECTION  3 
THEORY 

This  section  describes  the  theoretical  foundation  for  the  phased  array  lens 
simulator.  Although  the  theory  is  directly  applicable  to  active  lenses  composed  of 
arrays  of  metallic  radiators  (dipoles,  turnstiles,  etc.),  it  can,  with  straightforward 
modification,  be  extended  to  apply  to  other  lens  types  such  as  those  composed  of 
microstrip  arrays.  The  theoretical  development  is  not  complete.  Most  missing  de¬ 
tails,  however,  are  relatively  minor  and  are  expected  to  be  completed  and  implemen¬ 
ted  by  the  end  of  the  current  effort. 

A  theoretical  overview  is  presented  in  Section  3.1.  This  is  followed  (Section 
3.2)  by  a  detailed  discussion  of  the  application  and  generalization  of  infinite 
array  theory  and  the  moment  method  in  arriving  at  a  first  approximation  to  the  array 
element  currents.  The  general  theory  is  described  in  Sections  3.2.1  and  3.2.2  as 
it  applies  to  lenses  restricted  to  rectangular  array  lattices,  single-port  array 
elements,  planar  radiating  elements  (e.g.,  no  feedline  scattering),  and  a  special¬ 
ized  module  model.  The  theoretical  extensions  that  remove  these  restrictions,  model¬ 
ing  of  the  space  feed  illuminating  the  lens,  and  techniques  for  increasing  computa¬ 
tional  efficiency  are  discussed  in  Sections  3.2.3  through  3.2.8. 

The  theory  for  obtaining  corrections  to  the  infinite  array  currents  is  discussed 
in  Section  3.3.  A  higher  approximation  method  is  developed  that  employs  the  infi¬ 
nite  array  currents  as  the  initial  solution  in  an  iterative  procedure  that  accounts 
for  array  edges,  module  failures,  support  structure  scattering,  etc. 

The  theory  for  efficiently  computing  the  fields  radiated  by  the  array  element 
currents  is  contained  in  Section  3. A. 

3.1  OVERVIEW  / 

The  lens  array  simulator  solves  a  steady  state  (e^u't  temporal  variation  assumed) 
problem  by  a  successive  approximation  method.  This  method  is  particularly  appli¬ 
cable  to  lenses  composed  of  many  radiating  elements  (perhaps  hundreds  of  thousands) 
in  each  of  two  parallel  arrays.  The  first  approximation  is  based  on  infinite  array 
theory.  This  provides  approximations  to  such  parameters  as  radiating  element  cur¬ 
rents,  active  impedances,  coupling  between  arrays,  and  radiation  patterns  (by  trun¬ 
cating  the  infinite  array  after  the  currents  are  found).  The  resulting  radiating 
element  currents  are  used  in  obtaining  successively  better  approximations  of  these 
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parameters.  In  this  manner,  the  nominally  periodic  nature  of  each  of  several  sec¬ 
tions  of  a  lens  can  be  used  to  advantage  in  simplifying  the  overall  analysis.  The 
impact  of  such  effects  as  "feedback"  radiation  coupling  between  arrays  and  active 
impedance  variation  with  scan  can  be  assessed  with  infinite  array  theory  alone  since 
most  of  the  elements  are  expected  to  appear  as  if  they  reside  in  infinite  arrays. 

The  first  approximation  thus  is  of  use  in  lens  design;  e.g.,  for  optimizing  radiat¬ 
ing  element  type  (dipole,  folded  dipole,  etc.)  and  ground  screen  design  within  such 
constraints  as  module  load  impedance  limitations,  overall  lens  weight,  and  cost. 

3.1.1  First  Approximation 

Infinite  array  theory  is  applicable  to  large  numbers  of  identical  elements 
arranged  periodically.  A  candidate  lens  is  expected  to  be  largely  periodic;  however, 
due  primarily  to  manufacturing  and  deployment  constaints,  this  periodicity  is  ex¬ 
pected  to  differ  between  sections  of  the  lens.  In  the  first  approximation  method, 
since  most  sections  will  be  large,  infinite  array  theory  is  applied  separately  over 
each  section  —  called  a  "cell."  The  first  approximation  current  solutions  usually 
are  highly  accurate  within  each  cell;  they  may  be  less  accurate  near  the  cell  bounda¬ 
ries.  Higher  approximations  for  the  currents  on  these  elements  will  provide  corre¬ 
sponding  corrections  to  the  radiated  field  (Section  3.1.2).  Other  anomolies  in  peri¬ 
odicity  within  a  cell  due,  for  example,  to  element  module  failures,  also  will  be  re¬ 
solved  by  higher  approximations.  The  number  of  elements  requiring  analysis  beyond 
the  first  approximation  is  expected  to  be  relatively  few  (perhaps  hundreds  rather 
than  hundreds  of  thousands) . 

Each  cell  is  defined  by  a  periodic  lattice,  progressive  module  phasing,  uni¬ 
form  module  gains,  uniform  input  and  output  impedances,  and  identical  element  orien¬ 
tations.  Subdivision  into  cells  is  not  determined  solely  by  periodicity,  but  also 
by  relative  uniformity  of  module  characteristics  and  settings,  as  well  as  other  con¬ 
siderations  discussed  below. 

The  illumination  across  a  cell  during  radar  transmit  will  not  be  that  of  a 
uniform  plane  wave.  Application  of  infinite  array  theory,  therefore,  requires  plane 
wave  decomposition  of  the  illuminating  incident  field.  Each  plane  wave  excitation 
is  analyzed  and  the  resulting  currents  and  patterns  superimposed.  Since  a  component 
plane  wave  excitation  may  have  a  phase  distribution  across  the  lens  that  differs 
from  the  progressive  phase  setting  of  the  modules,  and  since  the  arrays  maysignifi- 
cantly  couple  via  RF  penetration  through  the  ground  screen  as  well  as  through  the 
modules,  "standard"  infinite  array  theory  alone  is  not  adequate.  (The  reference  to 


standard  infinite  array  theory  here  refers  to  an  analysis  of  a  periodic  array  of 
ident ical  structures  (elements)  —  either  plane  wave  or  lumped  circuit  source  ex¬ 
cited  —  where  the  latter  are  progressively  phased  along  the  array  plane.  Array 
elements  generally  are  not  identical  for  the  lenses  under  consideration  due  to  pro¬ 
gressively  phased  nodules.)  The  simulator  therefore  includes  an  extension  to  stand¬ 
ard  infinite  array  theory  whereby  "feedback"  modes  are  identified.  Each  feedback 
mode  assumes  element  currents  related  by  (3-1)  and  (3-2);  for  simplicity,  rectangular 
lattices  are  assumed. 

-ikmd  (s  +U' t  )  -iknd  (s  +ui  ) 
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Equation  (3-1)  applies  for  the  illuminated  array  (array  a);  (3-2)  applies  for  the 
nonilluminated  array  (array  b) 
where  (see  Figures  1-2  and  1-4) 

m  and  n  are  the  array  element  identifying  integers  (positive  and  negative 
v/here  m  =  n  =  0  refers  to  the  "reference"  element) 

d^  and  d^  are  the  x  and  y  interelement  spacings 

■  <x  and  -i  are  the  module  x  and  y  directed  progressive  phase  "directional 
cosine"  components 

s  and  s  are  the  incident  plane  wave  x  and  y  directional  cosines 
x  y  r 

u  is  the  feedback  mode  integer  (u  =  0,  1,  2,  ...) 

k  is  the  propagation  constant 

tX  b 

The  total  currents,  1^  and  I  ,  ate  obtained  by  superimposing  all  feedback  mode 
currents ;  i . e . , 
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The  feedback  mode  reference  element  currents  can  be  determined  iteratively, 
under  certain  conditions,  by  solving  first  for  the  zero1"'1  mode  currents  (u  =  0)  and 
so  forth  (Section  3.2.1).  Each  higher  order  mode  current  is  dependent  on  previous 
mode  computat ions ,  but  each  involves  essentially  standard  infinite  array  theory. 


That  (3-1)  is  a  reasonable  representation  can  be  argued  as  follows.  The 
array  a  current  phasing  is  influenced  by  the  phasing  of  the  incident  field.  The 
modules  impart  an  additional  phase  prior  to  inducing  array  b  currents.  The  result¬ 
ing  array  b  currents  radiate  through  the  ground  screen  and  induce  an  "additional" 
component  of  array  a  currents.  The  latter,  phased  by  the  modules,  induces  a  corres¬ 
ponding  component  of  array  b  currents  and  the  process  continues. 

The  standard  infinite  array  analysis  employed  is  based  on  a  plane  wave  expan¬ 

sion  formulation  for  the  field  from  an  infinite  array  of  current  dipoles  developed 
by  Munk  and  Burrell.  ‘  This  formulation  is  incorporated  into  a  "pulse  expansion 
function,  pulse  weighting  function"  moment  method  solution  for  the  element  currents. 
Radiating  elements  that  can  be  well  represented  as  identical  arbitrary  collections 
of  bent  wires  of  differing  radii  thus  can  be  analyzed  by  the  simulator.  The  major 

portion  of  the  computation  time  is  used  in  determining  a  corresponding  generalized 

impedance  matrix  [Z].'  This  matrix  relates  pulse  expansion  functions  in  typical  mo¬ 
ment  method  fashion  with  the  understanding  that  each  expansion  function  now  is  a 
progressively  phased  infinite  array  of  current  "pulses."  It  is  necessary  to  solve 
for  the  reference  element  current  only  since  the  unifoi.m  amplitude,  progressive 
phasing  of  the  element  excitation  implies  the  same  relationship  for  the  element 
currents.  The  matrix  equation  that  must  be  solved  is 

[  Z  ]  I  =  V  (3-5) 

where  the  coLumn  vectors  I  and  V  correspond  to  the  reference  element  current  distri¬ 
bution  and  excitation,  respectively.  The  number  of  simultaneous  equations  expressed 
by  (3-5)  is  limited  to  the  number  of  pulses  on  a  single  element;  thus,  "matrix 
solution"  consumes  little  computational  time.  By  far  the  most  computational  time 
is  used  in  computing  the  matrix  elements  of  [ Z ] .  Each  matrix  element  computation 
involves  a  truncated  doubly-inf inite  series  which  is  slowly  convergent;  also,  [Z]  is 
not  completely  independent  of  excitation..  The  matrix  must  be  recomputed,  for  exam¬ 
ple,  for  a  change  in  scan  angle. 

First  approximation  radiation  patterns  arc  computed  by  superimposing  the  ra¬ 
diation  from  each  feedback  mode  current  to  obtain  "total"  cell  patterns  and  then 
superimposing  these  cell  patterns.  Since  each  cellular  feedback  mode  current  is  a 
uniform  amplitude,  progressively  phased  distribution  —  as  indicated  by  (3-2)  — 
each  corresponding  radiation  pattern  is  computed  easily  by  pattern  multiplication. 
The  simulator  employs  a  closed-form  expression  for  the  array  factor  which  applies  to 
a  general  rectilinear  lattice  and  corresponding  cell  boundary. 


Lens  warp  tug  is  modeled  in  the  simulator  by  tilting  cells  to  best  fit  the 
curvature.  Since  each  cell  must  bo  planar,  the  accuracy  with  which  a  warp  is  to  be 
modeled  will  impact  on  the  choice  of  cell  size.  Greater  accuracy  usually  will  re¬ 
sult  from  larger  numbers  of  smaller  cells. 

3.1.2  Higher  Approximations 

First  approximation  currents  on  array  elements  near  cell  boundaries  and  dis¬ 
continuous  periodicities  (e.g. ,  nodule  failure)  within  a  cell  are  likely  to  be  in 
error.  Corrections  arc  provided  by  higher  approximations;  thus,  the  higher  approxi¬ 
mation  method  is  applied  only  in  the  neighborhood  of  cell  boundaries  and  other  dis¬ 
continuous  periodicities.  The  higher  approximation  method  can  include  the  effects 
of  scattering  from  support  booms  and  struts,  ground  screen  edges  and  creases,  etc., 
in  a  straightforward  manner. 

The  higher  approximation  method  makes  use  of  first  approximation  feedback 
fields  penetrating  the  ground  screen  and  the  fields  radiated  by  first  approximation 
currents  on  neighboring  elements.  These  fields,  in  addition  to  the  illuminating  in¬ 
cident  field,  form  a  known  excitation  from  which  the  second  approximation  current  on 

(?) a  (?)b 

an  array  a  element,  I  ,  and  the  current  on  its  array  b  counterpart,  1  ,  are 

computed.  This  computation  employs  a  "free  space"  thin  wire  moment  method  in  a 

straightforward  manner.  Since  the  radiation  coupling  through  the  ground  screen  is 

approximated  by  the  first  approximation  array  h  to  array  a  feedback  fields,  the  only 

coupling  between  the  arrays  that  must  be  solved  "simultaneously''  is  the  module  feed 

line  coupling.  The  latter  is  incorporated  easily  in  the  moment  method.  The  number 

of  neighboring  elements,  L,  that  must  be  considered  in  this  computation  is  determined 

(2)q  (7\|j 

by  observing  the  convergence  of  I'  and  I  as  L  is  increased.  This  number  is 
not  expected  to  be  large  (L  <  40)  for  reasonable  interelement  spacings  (i  A/2)  and 
most  scan  angles.  The  order  of  the  moment  method  matrix  equation  that  must  be 
solved  —  one  Cor  each  array  (a  and  b)  —  is  limited  to  the  number  of  expansion 
functions  on  a  single  radiating  element. 

The  third  and  higher  approximations  are  obtained  in  a  similar  manner.  Each 
approximation  uses  the  previous  approximations’ s  currents  to  obtain  the  contribution 
to  the  exciting  fields  from  neighboring  elements.  The  feedback  field  remains  the 
same  as  that  computed  in  the  first  approximation  except  for  a  scalar  multiplicative 
adjustment.  The  latter  is  determined  as  the  ratio  of  most  recent  port  current  approx¬ 
imation  to  first  port  current  approximation  for  the  array  b  element  under  considera¬ 
tion. 
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The  corresponding  higher  approximation  radiated  field  corrections  are  ob¬ 
tained  by  summing  the  fields  radiated  by  the  relatively  few  array  elements  that 
undergo  higher  approximation  current  corrections. 
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3.2  FIRST  APPROXIMATION  CURRENTS 


A  combined  infinite  array/moment  method  technique  is  applied  in  obtaining  a 
first  approximation  of  the  array  element  currents  in  a  cell.  (A  cell  is  defined  in 
Section  3.1.1  as  a  section  of  a  lens.)  This  method  is  an  extension  of  that  developed 
by  Munk  and  Burrell.1  It  is  applicable  to  a  large  lens  section  (cell)  composed  of 
identical  arbitrarily  bent  wire  radiating  elements  in  each  of  two  parallel  arrays; 
each  array  may  have  different  elements.  It  also  accounts  for  amplitude  and  phase 
adjust  modules  interconnecting  the  arrays  and  an  imperfect  ground  screen  between  the 
arrays.  In  particular,  the  technique  permits  the  modules  to  be  progressively  phased 
although  the  module  amplitude  adjustments  must  be  uniform  throughout  the  cell. 

3.2.1  Port  Representation  and  Solution 

A  lens  formed  by  two  infinite  dipole  arrays,  a  and  b,  in  periodic  lattices  on 

either  side  of  a  finitely  conducting  infinite  ground  screen  is  shown  in  Figure  1-2. 

The  two  arrays  and  ground  screen  are  parallel.  A  coordinate  system  with  x  and  y 

“*0X 

axes  located  on  the  screen  also  is  shown.  An  exciting  plane  wave,  E  ,  is  assumed 
arriving  from  the  negative  z  side  (behind  the  screen)  and  directly  incident  on  array 
a  (Figure  1-4);  amplitude  and  phase  adjust  modules  connect  array  a  to  array  b.  Array 
a  faces  the  radar  transmitter  during  radar  transmit;  array  a  faces  the  target  during 
receive.  Array  a  thus  is  always  on  the  illuminated  side  and  array  b  is  always  on 
the  nonilluminated  side. 

The  following  analysis  of  this  lens  is  not  restricted  to  dipole  array  elements. 
Each  element  can  be  a  collection  of  thin  bent  wires;  e.g.,  folded  dipole  or  dipole 
with  parasitic  scatterers.  For  simplicity,  all  wires  comprising  an  element  are  as¬ 
sumed  to  lie  in  their  respective  array  plane.  This  restriction  is  removed  in  Sec¬ 
tion  3.2.3,  where  elements  inclined  with  respect  to  the  array  plane  and  feed  line 
scattering  are  considered. 

Array  a  elements  may  differ  from  array  b  elements;  however,  each  element  within 
an  array  is  assumed  identical  and  lying  in  a  rectangular  lattice.  A  generalization 
to  arbitrary  rectilinear  (triangular,  etc.)  lattices  is  described  in  Section  3.2.4. 


The  analysis  is  structured  around  a  port  representation.  Each  pair  of  termi¬ 
nals  entering  an  array  element  is  assumed  to  be  a  port;  thus,  unbalanced  mode  cur¬ 
rents  at  each  terminal  pair  are  neglected.  The  module-to  radiator  feed  lines  are 
considered  part  of  the  modules.  Multiports  and  unbalanced  mode  currents  are  consid¬ 
ered  in  Section  3.2.7. 


The  port  voltages  and  currents*  are  shown  in  Figure  1-4.  The  superscript 
(a  or  b)  indicates  the  array  and  the  subscript  the  element  location;  e.g.,  mn  denotes 
the  mtb  column  from  the  yz  plane  and  ntb  row  from  the  xz  plane.  The  m  =  0  and  n  =  0 
elements  lie  on  the  z  axis  and  are  referred  to  as  the  reference  elements. 


Array  a  and  array  b  port  currents  are  determined  by  requiring  that  they  satis¬ 
fy  the  port  boundary  conditions 


va  =  vaa  +  vab  +  vex 

mn  mn  mn  mn 


(3-6) 


mn 


mn 
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where 

Va  and  Vb  are  related  to  la  and  Ib  via  the  mntb  module  (including  feed 
mn  mn  mn  mn 

lines)  two-port  parameters 

cici  t  h 

is  the  array  a  mn  element  port  voltage  with  F.  =  0,  the  ground  screen 

assumed  perfectly  conducting,  and  the  array  a  ports  excited  with  the  ia 

mn 

as  ideal  current  sources  (Figure  3-1) 

bfo  cici 

V  is  the  array  b  counterpart  to  V 

mn  mn 

Vab  is  the  array  a  mntb  element  port  voltage  with  EeX  =  0,  all  array  a  ports 
open-circuited,  and  the  1^  as  ideal  current  sources  exciting  the  array  b 
ports  (Figure  3-2) 

CX  til 

V  is  the  array  a  mn  element  port  voltage  with  all  array  a  ports  open- 
circuited  and  the  ground  plane  assumed  perfectly  conducting  (Figure  3-3) 

In  arriving  at  (3-6)  and  (3-7),  the  finite  conductivity  of  the  ground  screen 
is  considered  significant  with  regard  to  feedback  from  array  b  to  array  a.  Severe 
pattern  degradation  may  result  if  the  module  amplification  is  comparable  to  the 
ground  screen  attenuation;  therefore,  it  is  important  to  preserve  this  effect  in 
the  model. 

The  solution  of  (3-6)  and  (3-7)  is  based  on  techniques  recently  developed  by 
Munk  and  Burrell. 1  Each  array  element  is  modeled  as  piecewise  linear  and  the  array 
is  viewed  as  a  collection  of  infinite  arrays,  each  one  corresponding  to  one  linear 
segment.  Since  the  excitation  is  assumed  to  be  a  plane  wave,  the  relationship  of 


jwt 


*Phasors  associated  with  e 


temporal  variation  are  assumed  throughout. 


r 


Figure  3-2.  Problem  Defining  Vab 

mn 
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the  current  on  the  mnC'1  element  to  the  reference  (OO^1)  element  current  is  known. 

It  therefore  is  necessary  to  solve  for  the  reference  element  current  only.  A  moment 
method  procedure  can  be  applied  where  each  expansion  function  is  an  infinite  array 
of  currents  corresponding  to  one  linear  segment.  The  field  from  each  "expansion 
array"  is  expressed  as  a  doubly  infinite  sum  of  plane  waves.  One  advantage  of  this 
plane  wave  expansion  is  that  penetration  through  imperfect  ground  screens  becomes 
straightforward . 

The  method  presented  here  is  an  extension  of  Munk  and  Burrell's  work  in  that 
interconnecting  modules  are  introduced  between  arrays  and,  more  importantly,  the 
modules  may  differ  from  element  to  element.  In  particular,  the  modules  will  intro¬ 
duce  a  linear  progressive  phase  and  thus  differ  in  a  well-defined  manner.  The  mod¬ 
ules  are  identical  for  the  plane  wave  excitation 


E*X  (?) 


EeX 


(0) 


-jks^x  -jks  y 
e  e  y 


-jks  z 
z 
e 


The  induced  port  currents  become 


(3-8) 


Ia 

mn 


- j  kmd  s 

XX 


-iknd  s 

y  y 


00 


(3-9) 


,  ,  -jkmd  s  -jknd  s 

b  b  xx  YV 

I  =  e  e  J  y 

mn  00 


where 


r  =  xx  +  yy  +  zz  is  a  field  point 


(3-10) 


x,  y,  z  are  the  rectangular  coordinate  system  unit  vectors 
k  is  the  propagation  constant 

-*ex 

s  ,  s  ,  s  are  the  directional  cosines  of  the  propagation  direction  of  E 
x  y  z 

d^  and  d^  are  the  x  and  y  interelement  spacings 
3  b 

1^  and  Iqq  are  the  reference  element  port  currents. 

The  justification  for  (3-9)  and  (3-10)  can  be  demonstrated  from  linearity  considera¬ 
tions;  a  similar  form  applies  to  the  induced  port  voltages.  Equations  (3-9)  and 
(3-10)  are  not  sufficiently  general  if  the  modules  are  progressively  phased.  If  the 
phase  imparted  by  the  mn^  module  with  respect  to  the  reference  module  is 
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-k  (mdx  «x  +  nd^  a  ) ,  the  feedback  phenomenon  due  to  the  imperfect  ground  screen 
suggests*  port  currents  of  the  forms 


T  ia 
“  00(u) 


-ikmd  (s  +ua  )  -jknd  (s  +ua 

J  XXX  J 


to  run  v 

y  y  y) 


(3-11) 


,  .  -jkmd  rs  +(u+l)a  ]  -jknd  [s  +(u+l)a  1 

Ib  -  £  I*  ,  e  ■  x  x  x  e  y  y  y 

mn  u_q  00(u) 


(3-12) 


where  u  is  referred  to  as  a  feedback  mode  number.  The  Ia  .  .  and  Ib„,  .  are  the  u 

00(u)  00(u) 

mode  coefficients  of  array  a  and  array  b  reference  element  port  currents  determined, 
as  described  below,  by  satisfying  (3-6)  and  (3-7)  with  (3-11)  and  (3-12). 

With  port  currents  expressed  by  (3-11)  and  (3-12), Vaa,  Vbb,  Vab  in  (3-6)  and 

mn  mn  mn 

(3-7)  become 


-jknd  (s  +urt  )  -jknd  (s  +ua  ) 
a  a  x  xx  y  vy 

+ua  ,s  +ua  I00(u)  e  e 

x  x  y  y 


(3-13) 


Ebb 
Zs  + 


Zs  +(u+l)a  ,s  +(u-i-l)1'1  ^OO(u)  G 
x  x’  y  y 

-jknd  [s  +(u+l)a  ] 

y  y  y 


- j  kmdx [ sx+ (u+1) ax ) 


(3-14) 


,  ,  -jkmd  [s  +(u+l)a  ] 

ab  b  J  x  x  x 

“s  +(u+l)a  ,s  +(u+l)a  00(u) 

x  x  y  y 

-jknd  [s  +(u+l)a  ] 

,  y  y  y 


(3-15) 


z  is  the  array  a  active  impedance  when  driven  by  port  currents  (3-9)  and 

X  *  v 

J  the  ground  screen  perfectly  conducting 

zbb  is  the  array  b  active  impedance  when  driven  by  port  currents  (3-10)  and 
s  •  s 

X  V 

the  ground  screen  perfectly  conducting 


*  The  modules  inpart  a  phase  on  the  array  a  currents.  The  resulting  array  b  currents 
radiate  through  the  ground  screen  and  induce  an  "additional"  component  of  array 
a  currents.  The  latter,  phased  by  the  modules,  results  in  a  corresponding  com- 
- - -it  of  array  b  currents  and  the  process  continues. 


4 


zg  is  the  active  mutual  impedance  from  array  b  ports  to  array  a  ports  with 

X  ^  array  b  excited  by  port  currents  (3-10)  and  array  a  open-circuited 
(feed  lines  removed) 

Specifically, 


s  ,s 
x’  y 


o  .  S 

x  y 


s  ,  s 
x  y 


perfect  ground  screen 


perfect  ground  screen 


(3-16) 


(3-17) 


(3-18) 


imperfect  ground  screen,  array  a  ports  open-circuited 


where  IqQ  and  IqQ  are  reference  element  port  currents  with  the  arrays  excited  by  port 
currents  related  by  (3-9)  and  (3-10).  A  moment  method  technique  for  computing 

z  ,  z  ,  and  z  based  on  Munk  and  Burrell's  work,  is  described  in 

s  ,s  s  ,s  s  ,s  , 

x  y  x  y  x’  y 

Section  3.2.2.  An  expression  for  the  short-circuit  array  a  port  currents 
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Tex  Tex  xx  y  y 
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mn  00 


(3-19) 


due  to  E  with  the  ground  screen  perfectly  conducting,  is  determined  in  Section  3. 2 . 2 . 
ex 

The  V  in  (3-6)  then  are  given  by 
mn 


aa  _ex  aa 

z  I  =  -  z 

s  ,  s  mn  s  ,  s 

x  y  x  y 


-jkmd  s  -jknd  s 
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^00  e 


(3-20) 


where  the  reference  direction  for  I  is  the  same  as  for  I  in  Figure  3-1.  Equa- 

mn  mn 

ex 

tion  (3-20)  can  be  derived  by  considering  the  in  Figure  3-3  as  ideal  voltage 

sources.  The  current  flowing  through  these  sources  then  is  the  superposition  of  cur- 

“*ex  ex  gx  -»ex 

rents  due  to  E  with  the  V  sources  shorted  and  the  V  sources  with  E  removed. 

mn  mn 

Since  the  open-circuit  port  currents  are  zero, 

tex  +  '/!*  /  z*a  =0  (3-21) 

mn  mn  s  ,s 
x  y 

and  (3-20)  follows. 
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The  mnth  port  voltages  and  currents  are  assumed  related  by  the  mnth  module 


(including  feed  lines)  impedance  parameters  according  to 


„a  aa  Ta  ab  Tb 

V  =  -  z  ,  I  -z  ,  I 
mn  mod  mn  mod  mn 


,  ,  -ikmd  a  -iknd  a  , ,  . 
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V  =  -  z  ,  e  e  J  I  -  z  ,  I 

mn  mod  mn  mod  mn 


(3-22) 


(3-23) 


For  simplicity,  the  module  progressive  phase  setting  appears  only  in  the  array  a  to 
array  b  mutual  impedance  parameter.  General  module  models  are  considered  in  Section 
3.2.6.  In  combination  with  (3-11)  and  (3-12),  (3-22)  and  (3-23)  become 

“  -jkmd  (s  +ua  )  -iknd  (s  +ua  ) 
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,  ,  00  ,  -ikmd  [s  +(u+l)a  ]  -jknd  [s  +(u+l)a  ] 

_  bb  y-  b  J  xl  x  xJ  J  y 1  y  yJ 

mod  00  (u) 

u=0 

Equations  (3-6)  and  (3-7),  in  combination  with  (3-13),  (3-14),  (3-15),  (3-20), 


(3-24),  and  (3-25)  become,  after  dividing  through  by  e 


-ikmd  s  -jknd  s 

°  v  v  J  it 


ing  terms. 
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(3-26) 
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z  *  z  ,  +  z 
a  mod  s  +ua  ,  s  +ua 
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mod  I00(u)  +  Zu+1  ^(Xu)]  “  °  <3“27) 
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ab  _  ab  ab 

Zu+1  Zmod  Zs  +(u+l)a  s  +(u+l)a 
x  x  y  y 

bb  _  bb  bb 

zu+l  Zmod  Zs  +(u+l)a  ,s  +(u+l)a 
x  x  y  y 

Since  the  arrays  are  assumed  infinite,  corresponding  to  any  practical  module 
imparted  phase  shifts,  -kmd^a^  -  knd  a  ,  the  phasing  of  the  modules  repeats  every 
rectangular  section  of  M  by  N  elements;  i.e.,  M, M',N,  and  N’  can  be  found  such  that 


kd  a  , 
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(3-28) 


(3-29) 


Equations  (3-26)  and  (3-27)  then  become 
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=  0  (3-31) 


These  equations  must  hold  for  all  m  and  n.  Premultiplication  of  each  by 

j2rrmvM'/M  j27mvN'/N  ,  .  .  „  ,  , 

e  e  ,  where  v  is  an  integer,  summation  over  one  progressive  phase 

period"  1  <_  m  £  M  and  1  £  n  _<  N,  and  interchange  of  summations  results  in 
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Now 
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unless  (v-u)M’/M  and  (v-u)N'/N  are  both  Integers;  similarly 

ej27rm(v-u-l)M'  /M  ej2Trn(v-u-l)N'/N  _  q  (3-35) 

n=l  m=l 

unless  (v-u-l)M'/M  and  (v-u-l)N'/N  are  both  integers.  Assume  that  only  lower  ordered 
feedback  modes  are  significant,  as  is  expected  to  be  the  case  if  the  ground  screen 
attenuation  sufficiently  exceeds  module  amplification.  Then,  if  M  or  N  is  suffi¬ 
ciently  large,  it  follows  that  the  terms  in  (3-32)  and  (3-33)  containing  the  left- 
hand  side  of  (3-34)  can  be  assumed  nonzero  only  for  u  =  v  and  the  terms  containing 
the  left-hand  side  of  (3-35)  can  be  assumed  nonzero  only  for  u  =  v  -  1. 

The  solutions  to  (3-32)  and  (3-33)  then  are  determined  iteratively  as  follows. 

From  (3-32),  with  v  =  0, 
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From  (3-33),  with  v  =  1, 


(3-36) 


(3-37) 


ba  _a  bb  bb 

Zmod  00(0)  Z1  00(0) 


(3-38) 


X00(0) 


Zmod  Ta 
bb  00(0) 
Z1 


(3-39) 
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From  (3-32),  with  v  =  1, 
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And  so  forth;  in  general. 
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for  u  =  1,  2 .  In  conjunction  with  (3-11)  and  (3-12),  equations  (3-37), 

(3-44),  and  (3-45)  solve  the  problem. 
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3.2.2  Plane-Wave  Expansion  Moment  Method 

The  array  a  and  b  active  impedances,  z 


ab 


aa  bb 

,  z  ,  array  b  to  array  a  mu- 

S  •  S  S  i  s 

x  y  x  y 


tual  impedance,  z““  ,  and  the  array  a  reference  element  short-circuit  port  current, 

Sx,Sy 

GX 

Iqq,  are  considered  in  this  section.  The  development  is  based  upon  the  work  of  Munk 


and  Burrell. 


Consider  first  either  array  a  or  array  b  with  the  ground  screen  absent.  Let 
the  array  elements  be  composed  of  thin  wires  lying  in  the  plane  of  the  array. (The 
case  of  wires  with  components  normal  to  the  array  plane  can  be  handled  by  the  proce¬ 
dure  developed  below  if  some  formulas  are  replaced  with  those  of  Section  3.2.3). 

Let  the  array  lie  in  the  xy  plane.  A  port  voltage  or  current  source  array  excita¬ 
tion  that  is  uniform,  except  for  the  progressive  phase  -jk  (md^  s^  +  nd^  s  )  imparted 


at  the  mn^  element  with  respect  to  the  reference  element  (m  =  n 
the  current 


0) ,  results  in 


I  (r'  +  =  Wr,)  6 

mn  mn  UU 


-jkmd  s  -i  knd  s 
x  x  ^  y  y 


(3-46) 


th  ->■_ 

along  the  mn  element  wire  axis,  where  r  is  a  point  along  the  reference  element 

wire  axis  (path  C' )  and  r'  is  the  vector  from  the  reference  element  feed  to  the 
,  mn 

mn  element  feed.  This  form  for  I  (r  +r  )  also  will  arise  if  the  array  is  ex- 

mn  mn  1 

cited  by  a  plane  wave  with  directional  cosines  s  ,  s  ,  and  s  for  the  propagation 
direction.  Under  this  excitation,  once  Igg(r  )  is  determined,  the  current  on  any 
array  element  is  given  by  (3-46). 


The  boundary  condition  from  which  ^gg(r')  is  determined  is 


ES  =  Eex 
tan  tan 


tan 


(3-47) 


along  the  surface  of  the  reference  element,  where  E  is  the  field  radiated  by  all 


I  (r  +  r'  ),  E 
mn  mn 


tex¬ 


ts  an  exciting  plane  wave  in  the  absence  of  the  array  (impressed 
field),  E  is  the  total  field,  and  the  subscript  "tan"  denotes  "component  tangential 


to  wire  surface."  Along  perfectly  conducting  wires,  E 


tan 


0  at  all  points  except 


at'  a  source  port,  where  E  is  the  source  field  of  an  ideal  port  voltage  source, 

tsn  . 

Vqq,  and  it  is  assumed  that  each  mn  1  element  is  excited  by  port  voltage 


-jkmd  s  -jknd  s 

V  =  Vnn  e  X  X  e  ” 
mn  00 


(3-48) 


3-19 


A  moment  method  solution  of  the  thin  wire  approximation  of  (3-47)  results  in 
an  approximation  of  I00(r').2  In  this  procedure,  the  wire  axis  of  the  reference  ele¬ 
ment  is  approximated  by  an  Ng  straight  line  segment  fit (Figure  3-4).  The  reference 
element  may  be  composed  of  disjoint  wires  (as  with  parasitics)  and  also  may  have 
multiple  wire  junctions  (as  with  top-loaded  dipoles).  First,  Iqq (r ' )  is  expanded  as 


l'oo<;'>  •  z  b  b<;'> 
j  =  l 


(3-49) 


where  f.(r')  is  an  expansion  of  wire  current  constrained,  for  simplicity,  to  be  non- 
3  th 

zero  only  along  the  j  segment  of  length  ,  and  I.  is  an  unknown  coefficient  to  be 

J  J 

determined.  The  substitution  of  (3-49)  into  (3-47)  results  in 


/  ^  j  tan  tan  tan 


(3-50) 


where  iT^r.j)  is  the  field  radiated  by  the  inf inite  array  of  jth  segment  expansion 


_  ■jkmdxSx  'JkndySy 
functions,  f.(r')  e  e  y  .  With  respect  to  each  of  N  weighting  func- 

3  s  ” 

-F  — F 

tions,  wi(r),  the  weighted  average  of  (3-50)  along  a  path,  C,  on  the  reference  ele¬ 
ment  wire  surface  parallel  to  C' (Figure  3-4),  results  in* 


[Z ] X  =  vex  +  V 


(3-51) 
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where  the  i  elements  of  the  Ng  x  1  column  vectors  1,  V  ,  and  V  are,  respectively. 


=  /  rx.^(?)  <u 


(3-52) 


7.  =  f  E-w. 

1  JL 


4(r)  dJ, 


(3-53) 


and  the  ij  element  of  the  Ng  x  Ng  generalized  impedance  matrix  [Z]  is 


Etan(r.i)  di 


(3-54) 


*  For  convergence  purposes,  the  path  is  taken  a  wire  radius  in  front  of  the 
array  plane. 
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where  d l  is  an  increment  of  length  along  C.  The  reference  directions  for  the  paths 
C  and  C’  are  chosen  to  be  the  unit  vectors  and  i ^ ,  respectively,  where 
?..(?')  =  f^(r')jL  and  w^(r)  =  w^(r)i^  (Figure  3-4), 

In  determing  Z ,  the  field  Es  is  expanded  in  a  double  sum  of  plane  waves  as 
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for  (r-r)-z  <  0  where 
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r : 


th 


j  =  position  vector  to  the  center  point  of  the  j  segmert  along  C’ 
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free  space  wave  impedance 
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free  space  wavelength 


If  the  expression  under  the  radical  in  gz  is  negative,  gz  is  chosen  to  be 

negative  imaginary;  if  the  expression  is  positive,  gz  is  nonnegative  real.  This 

assures  that  all  plane  wave  modal  fields  travel  outward  from  the  array  plane  and 

nonpropagating  modes  attenuate.  A  plus  sign  subscript  on  h.,  g,  and  li!  corresponds 
_y  3  -  3 

to  field  points  r  in  front  of  the  array  plane;  i.e.,  (r  -  r!)*z  >  0.  A  minus  sign 

subscript  corresponds  to  field  points  behind  the  array  plane;  i.e.,  (r-rj  )*z  <  0* 

Equation  (3-55)  is  obtained  by  integrating,  with  respect  to  Z'  over  the 
segment,  the  expression  for  the  field  radiated  by  an  Infinite  array  of  point  current 
dipoles,  fj(r’)  d?,'  (Appendix  A).1  The  more  general  case,  where  ?j(r')  extends  over 
two  or  more  segments  bent  with  respect  to  each  other,  can  be  treated  in  a  similar. 
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although  more  complicated,  fashion.  Except  for  notations'!  difference,  (3-55)  fol¬ 
lows  the  development  of  Munk  and  Burrell.1 

■+  ■+  til 

If  w^(r)  is  restricted  to  extend  over  only  the  i  segment  along  C  of  length 

A^,  (3-54),  with  (3-55),  becomes 
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and  r^  =  position  vector  to  the  center  point  of  the  i  segment  along  C.  Only  the 
"+"  sign  is  indicated  for  g+  (hence,  also  for  ij>j+  and  ft  )  in  (3-57)  since  the  sim- 

plying  assumption  of  wires  lying  only  in  the  xy  plane  permits  choosing  r  as  well  as 
r^  with  no  z  components.  The  more  general  expression  for  Z_  is  derived  in 
Section  3.2.3. 


The  progressive  phasings,  s  and  s  ,  which,  for  a  particular  p  and  q,  cause 

x  y 

g  to  be  zero,  can  result  in  an  infinite  Z...  These  progressive  phasings  correspond 
z  ij 

to  "grating  lobe  singularities"  since  they  correspond  to  grating  lobes  just  entering 
the  visible  range;  i.e.,  nonattenuating  plane  waves  propagating  along  the  array 
plane.  This  is  one  possible  source  of  "blind  spots"  and  rapidly  varying  active  im¬ 
pedances  with  scan.  This  is  only  a  "potential"  condition  for  infinite  Z„  since, 
although  g  is  in  the  denominator  of  ft.  ,  in  some  situations  a  compensating  effect 

Z  J*T- 

occurs.  Such  is  the  case,  for  example,  in  E-plane  scans  with  periodic  infinite  ar¬ 
rays  of  y  directed  dipoles  not  widely  spaced  in  the  x  direction  since  ft^+ =  0  for 
p  =  0  and  q  such  that  gz  =  0  (grating  lobe  singularity).  The  presence  of  a  perfect 
ground  screen  potentially  can  inhibit  a  blind  spot  due  to  a  grating  lobe  singularity 
as  is  shown  below. 


If  the  array  is  a  distance  d  in  front  of  a  perfectly  conducting  infinite 
ground  screen  (positive  z  axis  directed  from  ground  screen  to  array  as  for  array  b), 
the  field  E  ,s  (r,j)  at  any  point  r  in  front  of  the  ground  screen  excited  by  the 
infinite  planar  array  corresponding  to  the  jt'1  expansion  function  and  its  image  is 
given  by 
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for  (r  -  rj)'z  <0.  The  first  term  in  (3-59)  is  the  ground  screen  absent  contribution 
given  by  (3-55).  The  second  term  is  contributed  by  the  perfect  ground  screen.  It 
is  important  to  note  that  £(  for  the  expansion  function  and  its  image  are  identi¬ 
cal  due  to  the  simplifying  constraint  that  all  wires  lie  parallel  to  the  xy  plane; 

i.e.,  d  is  independent  of  j.  Thus,  h.+  and  t|)!+  for  the  image  term  are  the  same  as 

3  3 

for  the  ground  screen  absent  term.  The  Ik  and  ij<j  do  not  occur  in  the  image  terms 
since  only  field  points  in  front  of  the  ground  screen  apply, 
til 

The  ij  element  of  the  generalized  impedance  matrix  for  the  array  in  front  of 

gs  ->s  gs  -><• 

a  perfect  ground  screen  [Z6  ]  is  given  by  (3-54)  with  E  ’ 6  replacing  E  ;  thus 
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As  with  (3-57)  only  the  "+"  signs  for  g+ ,  Ik  +  ,  andijj!+  appear.  A  rearrangement  of 
terms  in  (3-60)  results  in 
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It  is  apparent  from  (3-61)  that  Z^_.  is  bounded  at  a  grating  lobe  singularity  since 

sin  kdgz  bounded  as  g  -*0  (g  appears  in  the  denominator  of  h  ).  This  ground 

gz  Z  3 

screen  suppression  of  a  blind  spot  at  a  grating  lobe  singularity  assumes  the  radia¬ 
ting  elements  lie  in  the  array  plane.  If,  for  example,  the  dipole  arms  are  inclined 
with  respect  to  the  array  plane  or  feed  line  scattering  is  present,  the  ground  screen 
may  enhance  the  blind  spot  effect  at  a  grating  lobe  singularity  above  what  it  would 
be  if  the  ground  screen  is  absent;  thus,  consideration  of  out-of-plane  radiating  ele¬ 
ments  (Section  3.2.3)  is  important. 


With  a  ground  screen  present,  (3-51)  becomes 

[2gs]i  =  vex,gs  +  V 


(3-62) 


— cx  cs 

where  the  element  of  V  ,a  are  given  by  (3-52)  with  E  expanded  to  include  its 

th  — 6X  gs 

reflection  from  the  perfect  ground  screen;  thus,  the  i  element  of  V  ,e>  is  given  by 
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(3-63) 


The  shape  of  the  expansion  and  weighting  functions  ?.(r')  and  w.(r)  affect 

J  ^  -> 

only  1 1/'  and  it.  in  accordance  with  (3-56)  and  (3-58).  If  f.(r')  and  w.(r)  are  unit 
J  j±  i±  j  i 

"pulses,"  then 
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The  — x~— -  forms  in  (3-64)  and  (3-65)  exist  for  all  values  of  x;  however,  care  is 
x 

required  in  their  computation  when  x  is  near  zero. 

The  array  b  active  impedance,  z^  then  can  be  determined  by  solving  (3-62) 

x’  y 

GX  RS 

for  the  source  segment,  i  =  ig,  current  coefficient,  1^  ,  with  ,e  =0  for  all 


Then 


i  and  V.  =  0  for  all  1  4  i  • 
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The  array  a  active  impedance,  z3a  ,  can  be  obtained  in  the  same  manner. 

x*  y 

6X 

The  array  a  reference  element  short-circuit  current,  Iqq>  can  be  determined 

—  pg 

by  solving  (3-62)  for  I  .  Now,  however,  =  0  for  all  i  and  the  elements  of  V  ’ 


are  given  by  (3-63)  with  VaX  given  by  (3-52). 
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The  array  b  to  array  a  mutual  coupling,  z  ,  through  the  imperfect  ground 

Sx,Sy 

screen  can  be  determined  from  (3-51)  where  [Z]  is  the  array  a  generalized  impedance 

"~GX  -  ~^6X 

matrix,  =  0  for  all  i,  and  the  elements  of  V  are  given  by  (3-52)  with  E  re¬ 
placed  by  the  field  EC  transmitted  through  the  ground  screen  with  array  b  excited  by 
current  sources  given  by  (3-10).  (A  means  for  determining  E  is  discussed  below.) 


Thus 
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where  the  it'1  element  of  Vt  is 
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and  the  "a"  superscript  stands  for  "array  a."  The  i  element,  ,  in  the  solution 

s 

vector,  Ia  of  (3-67),  is  the  short-circuit  array  a  reference  element  port  current. 
The  reference  element  open-circuit  voltage  is 
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as  is  justified  by  the  argument  following  (3-20);  thus, 
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The  field  E  is  obtained  readily  from  knowledge  of  the  plane  wave  transmission 
coefficients  of  the  ground  screen  and  the  plane  wave  decomposition  (3-55).  The  array 
b  current  is  determined  with  array  b  excited  with  port  current  sources  (3-10).  The 
array  b  reference  element  current  is  obtained  from  the  solution  to  (3-47)  with 
V?X  =  0  for  all  i  and 
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The  solution  to  (3-51)  ,  denoted  X  ,  gives  rise  to  the  field 
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(3-72) 


incident  on  the  ground  screen  where  the  subscripted  negative  signs  indicate  that 
(3-72)  applies  for  (r-r! )-z<0  (which  is  appropriate  since  the  ground  screen  is  be¬ 
hind  array  b).  Since  E^  is  a  superposition  of  plane  waves,  E  is  readily  ob¬ 
tained  by  applying  the  proper  transmission  coefficients  to  each  component  plane  wave. 
Each  plane  wave  in  (3-72)  propagates  in  the  g_  direction;  thus,  g  and  z  (the  normal 
to  ground  screen)  form  the  corresponding  "plane  of  incidence."  A  unit  vector  normal 
to  the  incidence  plane  is 


g  x  z 


(3-73) 


and  a  unit  vector  parallel  to  the  incidence  plane  and  normal  to  g  is 
nn  =  g_  x  nj_ 


(3-74) 


Any  plane  wave  travelling  in  the  g_  direction  can  be  decomposed  into  two  plane  waves, 
one  having  the  E  field  polarized  along  nj^  and  the  other  having  the  E  field  polarized 
along  ii||  .  If  T^  and  T  |j  are  the  corresponding  transmission  coefficients  and  g_  is 

the  transmitted  wave  direction  (g  =  -z  for  a  good  conducting  screen J), 


3-27 


then 
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where  t|) .  is  with  g  replacing  g  . 

The  substitution  of  (3-75)  into  (3-68)  results  in 


'i  ~  X)  l]  2d^d  ^2  ^2 

j-i  p--“  q'-00 


-jk(r-?’)gt  r  ^ 

e  Lai-nL)Ti(nrhj_) 


(3-76) 


+  (L-ny  )T  ||  (n  || 
cib 

and  z  is  determined  from  (3-67)  and  (3-70)  as  previously  discussed.  In  (3-76) 

x’  y 

i|>^_  Is  i|>j_  with  replacing  g  . 

3.2.3  Inclined  Radiators  and  Feed  Line  Scattering 


The  development  in  Sections  3.2.1  and  3.2.2  is  limited  to  array  radiating 
elements  oriented  in  the  plane  of  the  array.  This  method  can  be  generalized  to 
handle  scattering  from  element  feed  lines  and  radiation  and  scattering  from  elements 
having  arms  inclined  with  respect  to  the  array  plane  (Figure  3-5) . 


Figure  3-5.  Coaxial-Fed  Bent-Arm  Dipole 


Munk  and  Burrell  [1,  eq.(3)]  obtain  the  expression  for  the  field  from  an  in¬ 
finite  planar  array  of  point  current  dipoles.  They  indicate  that  this  expression  is 
valid  for  arbitrarily  oriented  dipoles;  the  dipoles  need  not  be  oriented  parallel  to 
the  array  plane.  Munk  and  Burrell  also  mention  that  a  moment  method  can  be  developed 
from  this  expression.  This  development  is  carried  through  for  the  special  case  of 
array  elements  composed  of  bent  wires  confined  to  the  array  plane  (Section  3.2.2). 
Expressions  based  on  [1,  eq.(3)]  that  are  similar  to  the  mutual  impedance  expression 
[1,  eq.(19)]  given  by  Munk  and  Burrell  are  developed;  however,  [1,  eq.(19)]  was  de¬ 
rived  for  the  restricted  case  that  mutually  coupled  antenna  segments  did  not  both 
pass  through  any  plane  parallel  to  the  array  plane.  The  generalization  of  the  mo¬ 
ment  method  to  handle  "out  of  plane"  elements  necessitates  more  complicated  expres¬ 
sions.  These  are  derived  here. 

The  field  dE  at  r  due  to  an  infinite  planar  rectangular  array  of  point  cur¬ 
rent  dipoles 
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is  given  by 


-  r  -r  -r . 

dE(r ,r  )  = 


2d  d 
x  y 


E  E 


,-jk(r-r')-g+  £ 


P^  - 


oo  n=  —co 


h+f  d!. 


(r-r')  •zs-O 


(3-78) 


where  the  coordinate  system  is  that  of  Figure  1-2,  r'  is  the  location  of  the  refer¬ 
ence  element,  m  *  n  =  0,  and  the  remaining  variables  are  as  defined  in  Section  3.2.2 
with  replaced  by  the  unit  vector  £'.  This  result  is  derived  in  Appendix  A  for 
the  more  general  case  of  rectilinear  (skewed  grid)  lattices.  The  modifications  nec¬ 
essary  to  extend  the  results  in  this  section,  as  well  as  those  in  the  previous  two 
sections,  to  general  rectilinear  lattices  is  discussed  in  Section  3. 2. A. 

Consider  now  a  rectangular  array  of  identical,  but  arbitraily  bent,  thin  wire 
radiator  elements.  Let  the  reference  element,  m  =  n  =  0,  be  segmented  in  ac¬ 
cordance  with  the  moment  method,  as  in  Section  3.2.2,  and  let  fj(fc')  be  the 
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reference  element  segment  expansion  function.  Then,  with  r'  =  r^  +  as 

defined  In  Figure  3-6, 
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Implied  in  (3-80)  is  a  j  "  segment  integration  with  respect  to  A! .  This  integration 
divides  as  indicated  in  (3-82)  when  r  defines  a  point  common  to  a  plane  that  simul¬ 
taneously  divides  the  jth  segment  and  is  parallel  to  the  array  plane.  The  l]  co- 

^  J 

ordinate  of  the  point  along  the  j  segment  that  lies  in  this  plane  is  defined  by 

equating  z  axis  projections  of  r  and  r^  +  £(  f  t  ,  yielding 
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For  £j-z>0,  i.e.,  the  direction  indicated  in  Figure  3-6,  integration  with 

respect  to  from  -hjl  to  the  value  of  i,!  defined  by  (3-83)  involves  plane  waves 

propagating  away  from  the  array  plane  on  the  positive  z  axis  side;  hence,  the  +  sub¬ 
script  is  indicated  for  g  and  ti  in  the  first  integral  of  (3-82).  For  the  inte- 
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gration  from  the  (3-83)  value  of  f,l  to  Aj /2 ,  the  properly  chosen  plane  waves  propa¬ 
gate  along  the  negative  z  axis  away  from  the  array  plane.  Hence,  the  second  inte- 

„  ■+ 

gral  in  (3-82)  indicates  g_  and  hj_  with  negative  subscripts. 

The  -e  and  +e  appearing  in  the  first  and  second  integrals  in  (3-82)  are  in¬ 
cluded  to  avoid  integrating  through  the  singularity  that  occurs  in  the  plane  of  the 
source  array.  This  singularity  is  a  consequence  of  the  plane  wave  representation-  of 
the  field.  Although  this  singularity  appears  integratable,  this  adjustment  is  ex¬ 
pected  to  significantly  improve  convergence  of  the  doubly  infinite  sum  in  (3-80). 
For  nonself-elements,  i  4  j ,  e  equal  to  the  wire  radius  is  probably  adequate.  For 
self-elements,  i  =  j,  e  may  have  to  be  smaller. 

If  £  «z  <  0  (i.e.,  j  is  directed  along  the  jC^  segment  opposite  to  that 

indicated  in  Figure  3-6),  the  sign  of  (3-82)  simply  changes  as  indicated  by  the  ± 
depending  upon 
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As  discussed  in  Section  3.2.2,  the  ijc  generalized  impedance  matrix  element, 

tZ]±J ,  is  determined  by  weighted  integration  of  the  jth  current  segment  field  along 

the  surface  of  the  ith  wire  segment.  Let  r  locate  the  center  point  of  the  i  seg- 

^  th 

ment  surface  path  as  indicated  in  Figure  3-6).  (For  i  4  j ,  the  i  segment  match¬ 
ing  path  may  be  taken  along  the  axis  of  the  i*"*1  segment.) 
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where  is  given  by  (3-80).  For  simplicity,  assume  that  the  z  coordinate  projec¬ 
tions  of  all  wire  segments  are  either  disjoint  (nonoverlapping)  or  identical.  Then, 
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If  pulse  expansion  and  weighting  functions  are  considered,  then  Q  becomes 
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>*■ 

-jk(r 


r^>^«.e  i  r+ -e‘lk£i-^/2  J 

i  j+  jk£' -g+  L  +  i+J 


(3-87) 
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e  J  '  r  -e  ip 


,] 


if  )-z  =  0;  £j*z  ?  0 


where 


h±  =  Ai± 


sin|"k(£  ‘ g+) A  /2I 

kai-g±)A./2 


V  =  A 


sin[k(£’-g+)Aj/2] 

j  r(£’-i±)A./2 


sin  A+ 


A  -  k(£1*g±-£j-g±  £1-£/£^£)a±/2 
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If  a  ground  screen  is  present  a  distance,  d,  behind  the  array,  image  fields 
must  be  considered  in  each  impedance  element  computation.  Let  double  primed  posi¬ 
tion  vectors  and  segment  path  lengths  refer  to  image  currents.  Then,  since  the 
coordinate  origin  lies  on  the  ground  screen. 


rV  =  x(x*rj)  +  y(y-r!)  -  z(zt!) 


(3-88) 


i"  =  x(x'Jl!)  +  y(yl')  -  z(z-r'.) 

J  J  y  J 


(3-89) 


and  the  ijth  element  of  the  generalized  impedance  matrix  for  an  infinite  array  over 
a  perfect  ground  screen  becomes 


[*«*] 


rz]<4  + 


ij  1  Jij  2d  d 

J  J  X  y 


CO  00 

E  E 

p=  —oo  q=  -co 


3gs 

pq 


(3-90) 


where 


-jk(r  -r'')-g  gs 

0gs  =  e  1  J  +  n  -h"  il)  tj) 

Sq  1  j+  i+  J+ 


*  J+  '  J  k(q-g+)A./2 

For  plane  wave  incidence  excitation,  the  generalized  voltage  vector  elements 
given  by  (3-52)  must  be  appropriately  modified  to  account  for  reflection. 

3.2.4  Skewed  Lattices 


Many  array  lattices  are  skewed  rather  than  rectangular.  The  infinite  array 
analysis  described  in  this  report  is  applicable,  with  only  minor  change,  to  general 


rectilinear  (skewed)  lattices.  Any  point  in  a  skewed  point  lattice  can  be  located 
by  two  integers,  m  and  n.  The  position  of  the  mnth  point  is  given  by 

r  =  md  e  +  nd  y  (3-91) 

mn  e  y 

where  e  and  y  are  the  lattice  unit  vectors  with  associated  interelement  spacings 
d£  and  d^,  respectively  (Figure  3-7). 


A 

y 


Figure  3-7.  General  Rectilinear  (Skewed)  Lattice 


The  standard  x,v,z  rectangular  and  r,  9,<)>  spherical  coordinate  systems  are 

centered  at  the  m  =  n  =  0  point  and  oriented  as  shown  in  Figure  3-7.  If  identical 

th 

radiators  are  located  at  the  lattice  points,  the  mn  radiator  excitation  phase  that 
will  launch  a  wave  in  the  r  direction  is  expressed  by 

r-?'  =  md  s  +  nd  s  (3-92) 

mn  e  e  y  y 

where  s^  =  r-e,  =  r-y  and  the  indicates  a  unit  vector.  This  phasing  also  can 

be  expressed  in  terms  of  the  more  standard  x  axis  directional  cosine,  s^  =  r.x,  by 
noting  that  (see  Figure  3-7) 

d^e  =  d^x  -  Ay  y  (3-93) 


and,  with  (3-91), 

— ►  ^ 

r'  =  md  x  +  (nd  -  mAy)y 
mn  x  y  J  J 


(3-94) 


Hence, 


r-r'  =  md  s  +  (nd  -  mAy)s  (3-95) 

mn  xx  y  y 

It  then  can  be  shown 5  that  the  skewed  lattice  infinite  array  generalized  impedance 
matrices,  [Z]  and  [ZgS],  are  computed  by  the  formulas  developed  in  Sections  3.2.2 
and  3.2.3  if  the  expressions  for  g+  and  gz  are  replaced  by 


g.  3  s  r 

X  XA  y  I 

sx  +  p  T~  "  q  dd 

X  x  yj 

+  y 

[•x  +  ’r] 

T  x  y 

sx  +  p  d  ‘  q  Td 
x  x  y 

■2  r 

x  Xl2 
sy  +  q  r 
yj 

(3-96) 

(3-97) 


The  feedback  mode  development  in  Sections  3.2.1  and  3.2.6  also  applies  to 
skewed  lattices  if  the  phase  imparted  by  the  mnt'1  module  with  respect  to  the  refer¬ 
ence  module  is  expressed  by  -k(md  a  +  nd  a  ) . 

e  e  y  y 

All  d  ,  s  ,  and  a  in  these  developments  are  simply  replaced  by  d  ,  s  ,  and 
XXX  p  e 

.  .  i  aa  ,  bb  ,  ab 

a  ,  and  the  active  impedances,  z  ,  and  z  ,  and  z 

e  s  •  s  s*s  s 

e’  y  e’  y  e  y 

determined  with  the  expressions  for  g+  and  gz  given  above  where  the  relationship  be¬ 
tween  s  and  s  is  found  by  equating  (3-92)  and  (3-95)  resulting  in 
ex 


e  e 

are  understood  to  be 


s  = 
e 


d  S  -  AyS 

x  x  y  y 


(3-98) 
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3.2.5  Space  Feed  Model 

The  field  illuminating  the  lens  will  not  be  planar  during  radar  transmit; 
however,  a  plane  wave  decomposition  of  the  illumination  across  the  lens  face  is 
always  possible. 6  Each  plane  wave  component  of  the  illumination  can  be  applied  in¬ 
dividually  and  the  resulting  lens  patterns  superimposed.  The  first  approximation 
method  is  being  designed  in  this  manner. 

Since,  in  the  first  approximation,  infinite  array  theory  is  applied  to  each 
cell  individually,  a  separate  plane  wave  decomposition  of  the  illumination  is  car¬ 
ried  out  for  each  cell.  Only  a  small  number  of  plane  wave  components  must  be  con¬ 
sidered  for  each  case.  The  maximum  cell  size  may  be  chosen  such  that  a  single  plane 
wave  fit  to  the  illumination  across  a  cell  is  adequate  (Appendix  B) . 

The  propagation  direction  of  the  plane  wave  is  along  a  line,  from  the 
phase  center  of  the  feed  to  the  cell  center.  The  field  components  of  the  plane  wave 
are  the  projection  of  the  illuminating  field  at  the  center  of  the  cell  onto  the  plane 
perpendicular  to  P.  A  single  plane  wave  approximation  to  the  illumination  across  a 
cell  may  not  be  possible  if  a  large  number  of  cells  is  to  be  avoided.  A  10-percent 
accuracy  in  plane  wave  approximation  to  the  illumination  across  a  cell  could  result 
in  about  200  cells  subdividing  the  lens  (Appendix  B) .  A  tighter  accuracy  will  in¬ 
crease  this  number. 

3.2.6  Generalized  Module  Model 

The  idealized  module  model  assumed  in  the  Section  3.2.1  development  is  gen¬ 
eralized  here.  The  modules  were  represented  in  Section  3.2.1  by  two-port,  open-circuit 
impedance  parameters.  These  parameters  were  assumed  uniform  from  element  to  element, 
within  a  cell,  except  for  the  array  a  to  array  b  "mutual"  impedance  parameters  which 
differed  between  elements  in  progressive  phase  only.  In  this  section,  the  module 
parameters  are  permitted  to  be  any  complex  function  of  the  progressive  phase  setting. 
The  only  restriction  is  that  modules  with  phase  settings  differing  by  a  multiple  of 
2tt  have  identical  parameter  values.  Open-circuit  impedance  parameters  do  not  al¬ 
ways  exist;  e.g.,  a  half-wavelength  section  of  transmission  line  has  no  such  repre¬ 
sentation.  The  following  development,  therefore,  is  in  terms  of  two-port  scattering 
parameters.  The  assumed  form  of  the  "feedback”  mode  expansion  (3-12)  is  modified 
slightly.  This  permits  a  convenient  column  vector  representation.  The  resulting 
formulas  are  shown  to  reduce  significantly  for  the  special  case  of  a  passive  lens 
with  simple  line  length  phase  shifters. 
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The  n  i  element  array  a  and  array  b  port  currents  are  assumed  to  have  the 


-ikmd  (s  +  uot  )  -iknd  (s  +  ua  ) 
a  Jxx  x  Jyy  y' 

i00 (u)  6  6 


b  -jkmdx(sx  +  uax)  -jkndy(sy  +  uuy) 

I00(u)  6  6 


In  column  vector  notation,  (3-99)  and  (3-100)  can  be  represented  as 


-ikmd  (s  +  ua  )  -jknd  (s  +  ua  ) 

J  x  x  x  J  y  y  y  - 
e  e  j  j  J  lr 


I3  1 


The  port  conditions  that  must  be  satisfied  are 


vaa  +  vab 

mn  mn 


(3-99) 


(3-100) 


(3-101) 


(3-102) 


Let  the  mn  module  (plus  feed  lines)  scattering  parameter  matrix  be  denoted  [S  ]. 

mn 

a  b 

If  the  scattering  parameter  normalizing  impedances  are  R  and  R  ,  then 


V  =  [R]  '  ([U]  +  [S  ])C 
mn  mn  mn 


(3-104) 


i  =  [R]“1/2((U]  -  [s  ])c 

mn  mn  mn 


(3-105) 


where 


Ra  0 


0  R 


and  Ca  and  are  the  "incident"  scattering  variables, 
mn  mn  ° 

Assume  that  C  max-  be  expanded  in  "feedback"  modes,  as  is  I  ,  according  to 
mn  r  mn  6 


E-jkmdx(sx  +  uax)  -jkndy(sy  +  uuy) 

cc 


(3-106) 


Then  (3-103),  (3-104),  and  (3-106)  combine  to  yield,  after  division  by 


-ik(ms  +  ns  ) , 
e  x  y  ’ 


E-lkmd  ua  -jknd  uaj  , 

.  x  x,  •>  i  [R31/2au]  + 


[R!  am  *  ts„i>c00(u)  -  U„n00(u)  -  v“ 


(3-107) 


Another  equation  in  and  obtained  from  (3-101),  (3-105)  and  (3-106) ,  is 


-ikmd  ua 
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-iknd  ua 

y  y 


(R11/2((U]  -  [Sroj)c00(o) 


I00(u) 


0  (3-108) 


where  0  is  the  null  column  vector. 

A  set  of  simultaneous  equations  is  formed  from  (3-107)  and  (3-108)  in  a  manner  simi¬ 
lar  to  that  in  forming  (3-32)  and  (3-33).  The  module  phase  settings,  a  and  a  , 

x  y 

are  assumed  to  satisfy  (3-28)  and  (3-29);  thus,  fS  ]  repeats  every  M  x  N  elements. 

mn 

Premultiplication  of  (3-107)  and  (3-108)  by  e^2  mv^  ^  e~* 2^nvN  /N^  wjjere  v  is  an 
integer,  summation  over  one  "progressive  phase  period,"  llm^M  and  1-n-N,  and  inter¬ 
change  of  summations  yields 


([Fvu]C00(u)  "  tFvu]  [zJX00(u) 


)  = 


y_ ,  .-ex 

voo 


v=0 


v=+l,+2. 


(3-109) 


J2  1  500(u)  -  tFvul500(u)>  "  5  V=0’±;L’±2’  ••• 


(3-110) 


where 


N  M 
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vu 


^2^2  ej2m(v  “  u)M'/M  ej2  n(v  "  U>K'/N  tR]i5([U]-[Smn]) 


(3-111) 


n=l  m=l 

[*"  1  =  yy  ej27rm(v  -  U)M'/M  ej2TO(v  -  U)N'/N  [R]-^((U]  -  [S  ])  (3-112) 

vu  mn 
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otherwise 
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The  simplifications  resulting  in  the  right-hand  sides  of  (3-109)  and  (3-113)  assume 
that  (v-u)M'/M  and  (v-u)N’ /N  are  never  integers  simultaneously.  This  condition 
occurs  whenever  (v-u)  is  an  integer  multiple  of  M  and  N  simultaneously.  This 
condition  is  avoided  by  limiting  (3-109)  and  (3-110)  to  only  the  lower  valued 
equations  and  by  limiting  the  infinite  summations  to  only  lower  values  of  |u[ 

This  is  permissible  if  only  lower-ordered  "feedback”  modes,  u,  are  significant. 

A  special  case  of  interest  is  a  passive  lens  with  modules  replaced  by  trans¬ 
mission  lines.  The  relative  lengths  of  line  determine  the  progressive  phasing.  The 
scattering  parameters  for  the  mnth  element  become 

[S^]  =  -j 27rmM ’/M  sJ2toHVN  ^0  ^  (3-114) 

a  b 

if  R  =  R  =  Rg,  where  R^  is  the  characteristic  impedance  of  the  transmission  lens. 
Equations  (3-111)  and  (3-112)  reduce  to 


[F  *  ] 
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v  -  u  *  0 

v  -  u  -  1  =  0  (3-115) 

otherwise 
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v  -  u  =  0 

v  -  u  -  1  =  0  (3-116) 

otherwise 
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Equations  (3-109)  and  (3-110)  then  become 


^0  500(v)  +  ^  [WlSo(v-l)  “  ^zv]^00(v) 


*00  v=° 


0  V  =  ±1,±2 , 
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Finally,  (3-117)  and  (3-118)  can  be  combined  to  form 


(3-117) 


(3-118) 
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(3-119) 
v=+l,±2,  .. 


Equation  (3-119)  can  be  solved  for  under  the  previously  stated  condition  that 

only  lower-ordered  feedback  modes  are  significant.  Then  (3-118)  can  be  solved  for 

I00  (v)  ‘ 

3.2.7  Multiport  Array  Elements 

The  development  throughout  Section  3.2  has,  for  simplicity,  considered  single¬ 
port  radiating  elements.  The  extension  to  multiport  elements  (e.g.,  crossed  dipoles 
for  circular  polarization  (turnstile  antennas)  with  quadrative  phase  shifters  in  the 
feed  lines)  is  straightforward  due  to  the  moment  method  framework  employed 
[2,  Chapter  6].  The  present  simulator  development  effort  is  not  expected  to  imple¬ 
ment  this  extension;  however,  the  important  case  of  turnstile  array  elements  will 
be  modeled  by  transferring  the  quadrative  phase  shifters  from  the  feed  line  to  the 
appropriate  radiating  arms.  Turnstile  antennas  for  use  in  an  SBR  lens  array  are 
likely  to  be  designed  in  just  this  manner,  where  the  radiating  arms  are  either  re- 
actively  loaded  or  shaped  to  impart  the  desired  quadrature  phase  relationship. 

A  loaded  wire  can  be  modeled  easily  within  the  moment  method  framework.  If 

lumped  impedance  loads,  Z^,  j'  =  1,2, . . . ,N  ,  are  located  at  segments  j’on  a  wire,  the 

J  s 

voltage  column  vector,  V  in  (3-51)  becomes 


V  =  V’  -  CZ  I  I 


(3-120) 
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The  elements  of  the  column  vector,  V  ,  are  the  values  of  applied  voltage  sources  at 
driven  segments  and  [Z^]  is  a  diagonal  matrix  with  Zj,  for  its  j't}l  diagonal  element. 
Equation  (3-51)  then  becomes 

([Z]  +  [ZL]U  =  V6X  +  V'  (3-121) 

Equation  (3-121)  may  be  solved  for  I  as  discussed  in  Section  3.2.2. 

A  module  to  feed  line  junction  (module  output)  is  generally  a  three-terminal 
junction  with  the  third  terminal  necessitated  by  the  ground  screen.  Nonideal  balun 
arrangements  at  these  module  outputs  cause  unbalanced  currents  to  flow  on  the  feed 
lines  if  asymmetries  exist.  The  asymmetries  may  be  electrical  (e . g.  ,  module  phasing 
for  off-broadside  beam  steering)  or  structural  (e.g..  Figure  3-8).  The  net  feed  line 
current,  1^  +  1^ ,  is  counterbalanced  by  current  I  on  the  ground  -creen;  i.e., 

I  =  -  (I,  +  I„).  This  net  current  constitutes  the  unbalanced  mode  current.  If  I 
g  1  2  g 

is  significant,  a  module  output  generally  requires  a  two-port  representation.  (This, 
in  turn,  leads  to  a  four-port  representation  for  a  module  with  two  outputs.)  The 
initial  version  of  the  simulator  will  be  capable  of  modeling  a  module  output  only 
as  a  single  port.  The  effects  of  unbalanced  mode  feed  line  currents  can  be  ac¬ 
counted  for,  despite  this  limitation,  if  it  can  be  assumed  that  there  is  negligible 
coupling  within  the  module  between  the  balanced  and  unbalanced  mode  currents  and 
also  between  the  unbalanced  mode  currents  at  the  array  a  and  array  b  module  outputs; 
i.e.,  only  the  balanced  mode  currents  are  assumed  affected  by  the  module  circuitry 
and  the  feed  line  unbalanced  mode  currents  arise  entirely  from  scattering.  In  Figure 
3-9,  the  modules  are  treated  as  two  ports  (one  port  on  each  array  side)  and  the  cor¬ 
responding  port  voltages  and  currents  are  transferred  to  the  radiator  excitation 
ports  (indicated  by  the  sources)  via  the  feed  line  transmission  line  equations.  The 
composite  module/feed  lines  port  parameters  become  those  discussed  in  Sections  3.2.1 
and  3.2.6.  These  port  parameters  apply  strictly  to  the  balanced  mode.  With  radia¬ 
tor  excitation  ports  thus  related,  each  feed  line  is  modeled  as  a  single  thin  wire. 
The  feed  line  scattering  then  is  accounted  for  by  the  moment  method  technique  dis¬ 
cussed  in  Section  3.2.3.  Balanced  mode  feed  line  current  radiation  can  be  neglected 
due  to  the  close  spacing  of  the  wires  in  a  feed  line. 
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Figure  3-8.  Asymmetrical  Feed  Line  Configuration 
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b.  Triangular  (Piecewise  Linear) 


Figure  3-10.  Expansion  Functions 
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3.2.8  Increased  Computation  Efficienc 


Evaluation  of  the  generalized  impedance  matrix  [Z]  (or  [ 2®s ] )  of  Section  3.2.2 
usually  requires  considerable  computer  run  time.  This  is  especially  true  for  arrays 
of  thin  wire  radiators  oriented  parallel  to  the  array  plane  since  the  doubly  infinite 
summation  in  (3-57)  then  is  slowly  convergent.  Each  term  in  the  summation  corre¬ 
sponds  to  a  plane  wave  field  averaged  over  a  segment  a  short  wire  radius,  a,  in  front 
of  the  array  plane  containing  the  current  filaments.  The  plane  waves  corresponding 
to  large  |p|  or  |q|  decay  exponentially  away  from  the  array  plane.  The  arguments  of 


these  exponentials  are  proportional  to 


v'K 


J  +  [’  JJ 


where  d  and  d 
x  y 


the  x  and  y  coordinate  interelement  spacings;  thus,  wires  that  are  thin  with  respect 


to  d^  and  d^  result  in  slow  convergence.  Examples  indicating  the  number  of  p  and  q 


terms  (plane  wave  modes)  typically  required  are  given  in  Section  5.  Several  possibi¬ 
lities  for  reducing  the  number  of  times  that  (3-57)  has  to  be  computed  are  being  ex¬ 
plored.  These  techniques,  discussed  below,  also  apply  to  the  ground  plane  case 
represented  by  (3-61). 


The  Toeplitz  nature  for  straight  wires  is  being  employed.  The  [Z]  for  a 
straight  wire  divided  into  segments  has  2Mg  -  1  diagonals.  If  the  segments  are  of 
equal  length,  all  elements  of  [Z]  common  to  a  diagonal  are  equal;  i.e.,  [Z]  is  Toeplit 
The  H  “  elements  of  [Zl  are,  therefore,  determined  from  only  2Ns  -  1  computations  of 
(3-57)  —  one  for  each  diagonal.  If  an  array  radiating  element  is  a  collection  of 
straight  wires,  this  property  applies  only  to  submatrices  centered  along  the  main  di¬ 
agonal  of  [Z].  The  order  of  each  submatrix  equals  the  number  of  segments  on  the  cor¬ 
responding  straight  wire. 


Except  for  broadside  phasing  (s^=  =  0),  [Z]  is  not  symmetric;  however , 

[Z] . .  differs  from  [Z]„  only  by  a  sign  in  an  exponential.  This  permits  a  reduction 
of  almost  50-percent  in  [Z]  computation  time. 


The  development  in  Section  3.2.2  employed  "pulses  for  expansion,  f  (#.’), 

and  weighting,  w  (>’),  functions;  these  are  indicated  in  Figure  3-10a  .  Each  pulse 

occupies  one  wire  segment.  The  results  in  Section  5  indicate  that  many  pulses  are 

needed  to  correctly  compute  the  input  reactance  of  an  infinite  array  of  near  resonant 

thin  dipoles  —  approximately  20  per  half  wavelength  for  a  =  0.01A  and  approximately 

0.8^  interelement  spacing.  These  results  assume  equal  segment  lengths,  Aj  .  If  tiny 

A  are  required  only  in  the  feed  region,  a  large  reduction  in  N  may  be  possible  by 
J  s 

modeling  with  unequal  A..  Results  with  unequal  A  will  be  available  shortly. 


b.  One  Wire  Chosen 

Terminated  at  Junction 


Addition  of  a  Second  Wire  by  d.  Addition  of  a  Third  Wire  by 

Overlapping  First  Wire  Overlapping  First  Wire 


Figure  3-11 •  Modeling  a 


Junction  with  Overlapping  Expansion  Functions 


Overlapping  expansion  and  weighting  functions,  such  as  piecewise  linear 
(Figure  3-10b)  or  piecewise  sinusoidal  (Figure  3-10c) ,  can  be  tried  in  place  of 
pulses.  The  overlapping  of  adjacent  functions  assures  continuity  along  the  wire, 
possibly  resulting  in  a  need  for  fewer  expansion  functions  to  approximate  the  current 
than  are  required  with  pulses. 

Since  each  overlapping  expansion  function  traverses  two  segments,  each  element 
of  the  corresponding  [Z]  is  computed  by  the  sum  of  four  terms,  each  similar  to  (3-57) 
but  differing  in  computation  of  and  .  Care  is  required  when  treating  multiple 
wire  junctions  to  maintain  current  continuity  at  the  junctions.  This  can  be  accom¬ 
plished  by  following  a  simple  rule  when  modeling  wires.  To  illustrate,  consider  the 
three-wire  junction  shown  in  Figure  3-lla.  The  modeling  begins  by  choosing  any  one 
of  the  wires  as  terminated  at  the  junction  (Figure  3-llb).  The  remaining  wires  then 
are  added  in  succession  such  that  each  overlaps  any  previously  placed  wire  (Figures 
3-llc  and  d).  This  overlapping  is  accomplished  by  aligning  the  junction  point  with 
the  peak  of  the  end  triangle  (or  sinusoid)  function  on  the  wire  being  added.  This 
method  is  described  in  Reference  8  and  presented  in  terms  of  "independent  loop  cur¬ 
rents"  in  Reference  9. 

The  expressions  for  the  pulse,  triangle,  and  sinusoid  functions  are  given  in 
(3-122),  (3-123),  and  (3-124) (Figure  3-10). 

For  pulses, 

j  1  2  -  A./2  <  St'  <  2  +  A./2 

f.(2')  =  3  3  3  3  (3-122) 

3  I  0  otherwise 

For  triangles, 

1(2'  -  2.  .)/A7  2.  ,  i  St’  <  l 

3-1  J  3-1  i 

Uj+i  -  2  '  ) / Aj  l.  -  2'  ±  2j+1  (3-123) 

0  otherwise 
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For  sinusoids,  where  k  = 


sin  k(f  -  i  )/sin  kAT 
J  1  J 


,  -  !■'  <  i. 
J  1  J 


fj(r> 


sin  k(A  -  f)/sin  kA+ 
-1  3 


n, 

J  j+1 


(3-124) 


otherwise 


3.3 


HIGHER  APPROXIMATION  CURRENTS 


The  first  approximation  method  solves  for  the  ar.-ay  a  and  array  b  currents 
within  a  cell  assuming  the  cell  is  infinite  in  extent,  all  array  a  radiating  ele¬ 
ments  are  identical,  all  array  b  radiating  elements  are  identical,  and  interconnec¬ 
ting  modules  between  arrays  are  permitted  to  differ  only  in  phase  settings.  The 
variation  in  phase  settings  between  modules  must  be  linearly  progressive.  Each  cell, 
however,  will  be  finite  in  extent  and  may  exhibit  other  discontinuities  in  periodicity 
due,  for  example,  to  module  failures  or  nearby  scatterers  such  as  lens,  edges,  booms, 
hinges,  struts,  and  other  structural  supports.  The  first  approximation  currents 
for  elements  near  these  discontinuities  may  be  significantly  in  error.  These  cur¬ 
rent  predictions  can  be  corrected  by  the  higher  approximation  method  described  in  the 
following  subsections.  Although  there  are  a  great  many  elements  (hundreds  or  more) 
that  may  reside  near  discontinuities,  the  higher  approximation  method  is  expected 
to  prove  reasonably  efficient.  This  is  because  (1)  the  method  makes  use  of  the 
known  first  approximation  currents  and  feedback  fields  penetrating  the  ground  screen 
and  (2)  a  moment  method  matrix  has  to  be  computed  and  inverted  only  once  for  each 
type  of  radiating  element  (dipole,  folded  dipole,  etc.) in  the  entire  lens.  Each 
approximation  uses  the  results  of  previous  approximations.  The  required  number  of 
approximations  depends  upon  the  rapidity  of  convergence  of  successive  approximations. 

3.3.1  Port  Representation  and  Solution 

Let  the  first  approximation  currents  on  all  lens  radiating  elements  be 

known.  The  method  for  obtaining  higher  approximation  currents  on  an  array  a  and 

array  b  element  pair  near  a  discontinuity  in  periodicity  is  presented  below  by  con- 

sidering  the  second  approximaton  and  extending  the  results  to  the  n  approximation. 

A  single-port  representation  of  a  radiating  element  is  assumed.  The  method  can  be 

expanded  to  include  multiport  models.  The  problem  being  modeled  by  the  second  approx 

imation  method  is  shown  in  Figure  3-12.  The  second  approximation  currents  on  a  lens 

array  a  and  array  b  element  pair  are  sought .  For  discussion  purposes,  these  elements 

are  referred  to  as  the  zero*"^1  elements  and  variables  corresponding  to  this  element 

pair  are  subscripted  zero.  Nearby  element  pair  currents  (subscripted  1,  2,...)  that 

t  h 

are  considered  in  computing  the  zero  element  pair  currents  may  reside  in  the  same 
cell  or  nearby  cells.  A  parenthetical  number  appearing  in  a  superscript  of  a  vari- 

(2)a 

able  denotes  the  level  of  approximation  of  the  variable;  e.g.,  is  the  second 

approximation  port  voltage  of  the  zero1"'"'  array  a  element,  and  I^^3  is  the  first 
approximation  oort  current  on  the  "second"  of  the  array  a  "neighboring"  elements. 
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.j.  m  .1  ..  .  I. 


The  incident  exciting  field  (space  feed  on  radar  transmit,  target  return  on  radar 

->-0X 

receive)  is  denoted  E  .  The  second  approximation  port  voltages  and  currents  for 
the  element  pair  under  consideration  (zero11*1  element  pair)  are  related  by 


v(2)a  =  a  (2)a  ex  y(l)ab  (l)aa' 
0  0  0  0  0 

(2)b  =  b  (2)b  (l)bb' 

0  0  0 


(3-125) 

(3-126) 


where 


a,  b. 
z  (2  ) 


V 


(l)aa' 

0 


is  the  impedance  of  an  isolated  array  a  (b)  element  with  the  ground 
screen  perfectly  conducting 

(V<1)bb)  is  the  array  a  (array  b)  2erotb  element  open-circuit  port 
voltage  due  to  the  nearby  array  a  (array  b)  first  approxi¬ 
mation  currents  with  the  ground  screen  perfectly  conducting 


„(l)ab 

0 


is  the  array  a  2erotb  element  open-circuit  port  voltage  due  to  the 
first  approximation  array  b  excited  field  penetrating  the  ground 
screen 


6X  th  ->ex 

Vq  is  the  array  a  zero  element  open-circuit  port  voltage  due  to  E  in 
the  presence  of  the  ground  screen 

(  2  )  s 

The  reference  directions  for  the  open-circuited  voltages  are  the  same  as  for  Vg 

and  The  za  and  zb  "isolated  element"  impedances  can  be  determined  with  feed 

line  scattering  accounted  for  as  suggested  by  the  Figure  3-9  model.  The  feedback 

voltage  Vg1  ab,  is  available  from  the  first  approximation  analysis. 

The  quantities  Vg^a,  Vg^b,  Ig^a  ,  and  Ig^b  are  related  by  module  two-port 

parameters;  e.g.,  the  impedance  parameter  representation  is 


,(2)a  aa  T(2)a  _  ab  T(2)b 
0  Z0  0  Z0  0 


.,(2)b  _  ba  (2)a  bb  _(2)b 
0  “  0  i0  “  0  0 

Equations  (3-125)  through  (3-128)  can  be  combined  to  form 


(3-127) 

(3-128) 


f  a  +  T(2)a  ab  (2)b  ex  (l)aa'_  v(l)ab 

U  z0  ’  0  0  i0  v0  V0  0 


(3-129) 
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m 


ba 

Z0 


(zb  +  zbb)  I(2)b  =  -  V(1)bb' 
(.z  +  z0  ;  J-Q  VQ 


(3-130) 


These  equations  explicitly  relate  second  approximation  port  currents  and  first  approx¬ 
imation  port  voltages. 

A  thin-wire,  pulse  expansion-impulse  weighting,  free-space  moment  method  code 
can  be  used  to  determine  z3,  zb,  and  V3X  and,  from  first  approximation  currents  and 

fields,  Vgb^aa,  Vgb^bb,  and  VQb^ab.10  Other  moment  method  codes  can  be  used  (e.g., 

see  Refs.  8  and  9);  however,  the  pulse  expansion  code  fits  well  with  the  pulse  ex¬ 
pansion  infinite  array  technique  employed  in  the  first  approximation  method  of 
Section  3.2.2. 

Moment  methods  calculate  short-circuit  currents  more  efficiently  than  open- 

.  ,  ,  ,,ex  ..(l)aa'  „(l)bb'  ,  ,.(l)ab  .  ,  „  th 

circuit  voltages;  therefore,  Vq  ,  ,  V  and  Vq  are  determined  from  zero 

element  array  a  short-circuit  port  currents  I^X,  ,  and  I^b^bb  and  zerotb  element 

array  b  short-circuit  port  current  I^b^ab  according  to 


..ex  a 

V0  -  -  z 

(3-131) 

(l)aa'_ 

v0 

a  T(l)aa' 

z  xo 

(3-132) 

V^bb'=  - 
0 

b  T(l)bb' 
z  xo 

(3-133) 

v(1)ab  =  - 
0 

a  T(l)ab 

z  xo 

(3-134) 

The  reference  directions  for  these  short-circuit  currents  are  the  same  as  those  for 
Ig2^a  and  I^2^b.  The  justification  for  (3-131)  through  (3-134)  follows  closely  the 
argument  following  (3-20). 

In  the  pulse  expansion  function  moment  method,10  an  isolated  element  is  mod¬ 
eled  as  Ng  connected  segments  (Figure  3-4).  Let  i  be  the  number  of  the  segment 
which  contains  the  element  port.  A  corresponding  isolated-element-above-a-ground- 
screen  "generalized  admittance  matrix"  can  be  determined. 0  Each  short-circuit  port  cur¬ 
rent  in  (3-131)  through  (3-134)  can  be  determined  from  row  ig  of  the  appropriate  gen¬ 
eralized  admittance  matrix  and  the  appropriate  "generalized  voltage  vector"2  accord¬ 
ing  to 
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Tex  'a  -ex 
I0  Yi 

s 


T  (l)aa' _  ~a  -(l)aa' 

0  i  0 

s 


T(l)bb’_  -b  =(l)bb' 

0  i  0 

s 


(l)ab  'a  - ( 1 ) ab 

0  i  v0 

s 


(3-135) 


(3-136) 


(3-137) 


(3-138) 


In  (3-135)  through  (3-138),  Y^  (Y^  )  Is  the  matrix  row  vector  taken  from  the  i  row 

s  s 

of  the  array  a  (array  b)  isolated-element-above-a-ground-plane  generalized  admittance 

matrix,  v^^aa  (Vgb^bb  )  is  the  array  a  (array  b)  zerotb-element-above-a-ground-screen 

generalized  voltage  vector  resulting  from  the  nearby  array  a  (array  b)  first  approxi¬ 
mation  currents,  is  the  array  a  zerotb-element-above-a-ground-screen  general¬ 

ized  voltage  vector  resulting  from  the  first  approximation  field  penetrating  the 

->£  _ 0  X  til 

ground  screen  (E  in  Section  3.2.2),  and  V  is  the  array  a  zero  _element-above-a- 

D-gx 

ground-screen  generalized  voltage  vector  resulting  from  E 
From  (3-129)  through  (3-138) , 


(z3  +  zaa)l(2)a  +  zab  I(2)b  =  za  Ya  (V6X  +  V(1)aa'+  V(1)ab) 

(  0  ;  0  0  0  i  '•  0  0  0  ' 

S 

ba  (2)a  b  zbb}  (2)b  =  b  -b  -(l)bb' 

0  0  u  0  '  0  i  0 

S 

Equations  (3-139)  and  (3-140)  can  be  solved  for  I^2^3  and  lQ^b.  1 


(3-139) 


(3-140) 


The  second  approxi¬ 


mation  zero  element  active  impedances  are  determined  from  (3-127)  and  (3-128)  as 

(2)a  .T(2)a  f  ,  ,,(2)b/T(2)b  ,  . 

/Iq  for  array  a  and  /Iq  for  array  b. 

Many  elements  are  likely  to  reside  near  discontinuities  in  periodicy.  This 
implies  a  large  number  of  (3-139)  and  (3-140)  computations  —  one  for  each  of  these 

**g 

elements.  The  admittance  row  vectors,  Y^  and  Y^  ,  are  fixed  for  all  (3-139)  and 

s  s 

(3-140)  computations  if  all  the  radiating  elements  throughout  the  lens  face  are  identi- 
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cal. 


This  also  applies  to  z3  and  z\  since  za(z^)  is  the  reciprocal  of  the  ig  ele¬ 
ment  of  Ya  (Y^  ).  The  major  effort,  then,  is  in  computing  V3X,  Vq^33  ,  ,  and 

s  s 

^ (l)ab 
0 

Equations  similar  to  (3-139)  and  (3-140)  can  be  found  for  the  array  a  and 
array  b  second  approximation  current  distributions.  These  current  distribution  equa¬ 
tions  employ  the  entire  isolated  element  generalized  admittance  matrices,  from  which 
~ci  ~b 

Y.  and  Y.  are  extracted,  with  significant  increase  in  computation  time.  The  simu- 
s  s 

lator  is  being  designed  assuming  the  first  approximation  current  distributions  ade¬ 
quately  describe  the  current  distributions  for  any  higher  approximation;  i.e. ,  a 
higher  approximation  current  is  assumed  to  differ  from  the  first  current  by  a  complex 
constant  scaling  factor. 

The  third  and  higher  port  current  approximations  are  obtained  in  a  manner 
similar  to  that  for  the  second  approximations.  The  zeroth  element,  nth  approxima¬ 
tion,  array  a  and  array  b  port  currents,  Iq3'*3  and  Iq”^,  are  obtained  from  the 
n-1  approximations  by  solving 


f  a  +  7aa’iT(n)a  ab  (n)b  a  ~a  ,^ex  -(n-l)aa'  -(n-l)ab 
U  0  M0  0  0  i  iV0  0  0 

S 


(3-141) 


ba 

Z0 


T (n ) a  ,  /  b  bb.  (n)b  _  b  ~b  -(n-l)bb' 
1Q  +  (z  +  z0  )I0  -  z  Y.  V0 


(3-142) 


s 

The  free  space  field  at  the  zeroC^  array  a  element  arising  from  the  n-1  approximation 
field  penetrating  the  ground  screen  results  in  V^3  This  field  is  difficult  to 

compute  exactly.  The  approximation  resulting  from  equating  V^3  and  is 

being  considered.  Another  possibility  is  "scaling"  according  to 


T  (n-l)b 

.-.(n-l)ab  _0 _ *(l)ab 

0  T(l)b  V0 

0 


(3-143) 


_  i  n(n-l)aa'  ,  r.(n-l)bb' 

The  column  vectors  and  Vq 

from  L  nearby  elements;  i.e.. 


involve  summations  of  the  fields 


L 

-(n-l)aa'  =  -(n-l)a2 

0  La  0 

2=1 


(3-144) 
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^n_1)bb’=  j2~vt1)hl  (3-iA5> 

1=1 

where  (V^11  b^bS-)  is  the  contribution  to  the  zerotb  element  array  a  (array  b) 

generalized  voltage  excitation  vector  from  the  n-1"*1  approximation  current  on  the 
th 

i  ‘  nearby  element.  These  computations  are  facilitated  by  assuming  the  first  approx¬ 
imation  distribution  of  current  on  all  nearby  elements  is  adequate  for  all  approxima¬ 
tions  and  scaling  in  accordance  with  the  higher-approximation  port  currents;  i.e.. 


=(n-l)aa' 

0 


(n-l)a 

± 

(l)a  0 
1S. 


(3-146) 


y(n-l)bb' 

0 


£ 

0  =  1 


I (n-l)b 

j(l  )'b 
Li 


V<x>w 

0 


(3-147) 


(n)a  (n)b.  .  ,  .  th 

1^  (I^  )  is  the  array  a  (array  b)  n  approximation  port  current  at  the 

nearby  element. 


i 


th 


3.3.2  Support  Structure  Scattering 

The  effects  of  hinges,  ground  screen  edges,  and  other  support  structures  on 
the  radiating  element  currents  can  be  accounted  for  naturally  within  the  higher  ap¬ 
proximations.  The  first  approximation  field  excites  scattering  currents  on  the 
support  structures.  These  currents  radiate  fields  that  enter  into  the  higher  approx¬ 
imations  by  the  addition  of  generalized  voltage  excitation  vectors,  VaS  and  VqS,  to 
(3-141)  and  (3-142).  These  equations  become 


a  aa.  (n)a  ab  (n)b 
+  z „  )I„  +  z„  I. 


/  a  . 

(z  +  z0  K0 


Jo  0 


ba  (n)a  ,  b  bb.  (n)b 
z0  l0  +  (z  +  z0  }I0 


za  Ya  (V-  +  vjn'1)aa  +  vjn-1)ab  +  Vas)  (3-148) 
s 

zb  ~b  ^(n-l)bb’  +  v<“-DbS)  (3-3.49) 

s 


The  infinite  array  theory  employed  in  the  first  approximation  method  naturally  pro¬ 
vides  a  plane  wave  spectra  representation  of  the  fields  exciting  the  support  struc¬ 
tures.  The  numerous  moment  method  and  ray  theory  computer  codes  available  then  are 
usually  directly  applicable  to  computing  the  support  structure  scattering. 
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3.4  RADIATION  PATTERN 

Radiation  patterns  are  determined  from  the  target  side  radiating  element  cur¬ 
rents  on  array  a  (array  b  with  radar  in  transmit  mode).  If  possible,  reciprocity  is 
emploved  during  radar  receive.  If  reciprocity  is  not  permissible,  as  is  likely  with 
active  modules,  the  radar  receive  pattern  has  to  be  analyzed  directly.  The  array 
b  radiating  elements  then  lie  on  the  feed  side  of  the  lens.  Far-field  radiation  for¬ 
mulas  then  apply  to  lens  sections  individually, provided  each  section  is  in  the  far 
field  of  the  receiver.  In  many  instances,  the  cells  defined  in  Section  3.1.1  may 
satisfy  this  constraint. 

The  first  approximation  currents  on  a  cell  have  associated  closed  form  radiation 
patterns.  These  are  employed  where  ever  possible.  The  details  in  computing  first 
approximation  cell  and  lens  far-field  patterns,  with  and  without  lens  structural 
warping,  are  discussed  below. 

3.4.1  Planar  Periodic  Array 

Radiation  from  one  side  of  a  cell  having  boundaries  coinciding  with  its 
periodicity  axis  (Section  3.2.4)  is  considered  here.  Cells  of  more  general  shape  are 
treated  in  Section  3.4.4.  All  radiating  elements  throughout  the  cell  are  identical, 
all  are  located  at  the  nodes  of  a  rectilinear  lattice  (Figure  3-7),  all  lie  a  height 
d  above  a  perfect  ground  screen,  and  all  can  be  modeled  as  electrically  thin  bent 
wires.  A  typical  cell  and  coordinate  system  are  shown  in  Figure  3-13.  The  position 
vector  r'  locates  a  point  on  a  radiating  element,  r"  locates  its  image  point,  and  r 
locates  a  field  point  P.  The  coordinate  system  xy  plane  coincides  with  the  infinite 
ground  screen.  The  far-field  vector  potential  for  this  array  is  given  by 


A(r)  -=  e 


-jkr 


4irr 


/ 


I'(r')  e 


jk- 


dr'  + 


fv 


(r”)  e 


jk- 


r  •  r 


dr" 


(3-150) 


where  I'(r')  is  the  array  current,  I"(r")  is  the  image  current,  and  k  is  the  wave- 
number.  The  field  point  direction  is  given  by 

r  =  sin  0  cos  <!>  x  +  sin  0  sin  y  +  cos  0  z  (3-151) 

where  a  indicates  a  unit  vector.  The  xyz  components  of  r  are  the  directional 


cosines  of  r. 


The  array  current  and  image  are  sums  of  the  N  individual  radiating  ele¬ 
ment  currents  and  their  images  IV. 

Equation  (3-150)  then  becomes 

(3-152) 


A(r) 


-jkr 


4nr 


±  r/i 


(r’) 


jk£- 


r . 
1 


dr’  + 


fVl 


(r") 


jkr-r 


dr" 


i=l 


where  r!^  locates  a  current  source  point  on  the  i  element  and  rV  locates  its  image 
point;  i.e.. 


1 


-V  •>  ^  ^ 

(rl  •  x)x  +  (r'  •  y)y  -  (r|  •  z)z 


as  indicated  in  Figure  3-14.  Let 


(3-153) 


r!  =  r  .  +  t' 

l  oi 


r  +  t" 
or 


(3-154) 

(3-155) 


where  r^  locates  the  ground  screen  point  corresponding  to  the  array  lattice  point 
(the  module  is  generally  located  at  rQi  ),  and  t'  and  t"  are  as  indicated  in  Figure 
3-14.  Note  that  t"  is  defined  from  t'  by 


t"  =  (t'  •  +  (t'  •  y)y  -  (t'  •  £)£ 


(3-156) 


Also  t'  is  the  same  for  corresponding  points  on  the  structurally  identical  elements. 
Equation  (3-152)  becomes 

(3-157) 


A(r')  =  e 
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4irr 


jkf'roi 


i=l 


/q<; 


’)  e 


jkr*t ' 


*h 


dr’  +/ I"(r")  e 


jkr  •  t 


Assume  that  the  current  distributions  on  all  elements  throughout  the  cell  are  identi¬ 
cal  within  a  complex  constant;  i.e., 

i:  =  K.r  ,  (3-i58) 

i  l  ref 

th 

where  1'  ,  is  the  reference  element  current  and  K.  is  the  ratio  of  i  element  to 
ref  i 

reference  element  port  currents. 
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Equation  (3-157)  becomes 


(3-159) 


A(r)  =  e 
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Eh  •  01 
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where  I"ej(r")  is  the  reference  element  image  current. 

The  reference  element  current  moment  can  be  approximated  by 


N 

s 


E 


Tref . 
1 


a.ji. 

i  j 


where  the  electrically  small  segment  lengths  A.  and  segment  unit 


vectors  &!  are  shown  in  Figure  3-4,  and  N  is  the  total  number  of  segments  on  the 
J  s 

wire  (note  that  the  segmentation  begins  a  half  segment  in  from  each  wire  end;  this 
better  models  the  zero  current  condition  at  these  points).  The  corresponding  image 


N 

s 

current  moment  becomes  -  /  I'  -  A.i." 

ref  j  j 

j=l  J 


where 


i”  =  (£j  '  *)x  +  Wy  y)y  -  (£j  •  z )i  (3-160) 

Equation  (3-159)  then  is  approximated  by 

A(r)  =  eF  (3-161) 

where  the  element  factor  e  and  array  factor  F  are  given  by 


-jkr 


4irr 


j  =  l 


1 '  ,  A. 
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(3-162) 


v — .  jkr-r  . 

*  ■  E  «.  •  01 

i=l 


(3-163) 


-•>  ->  ->  -> 

respectively.  The  t!  and  t.  in  (3-162)  correspond  to  t*  and  t"  at  the  center  points 
t  i  1  3 

of  the  j  segment  and  its  image,  respectively.  Equation  (3-161)  is  an  expression  of 
pattern  multiplication.  The  underlying  assumption  in  (3-161)  is  that  the  current 
distributions  on  any  two  array  elements  differ  by  only  a  complex  constant. 

The  6  and  if  components  of  the  far  radiated  E  and  H  fields  are  given  in 
terms  of  F  and  the  6  and  components  of  e  as 


E0  =  -j^eQF 

% =  -JknV 

"e  -  -  t 

%  =  "  n  Ko 

where  n  is  the  free  space  impedance. 


(3-164) 


3.4.2  Regular  Cells 

The  array  factor  F  assumes  a  closed  form  expression  for  each  feedback 
mode  of  the  first  approximation  currents  on  a  regularly  shaped  cell.  A  regularly 
shaped  cell  is  one  with  boundaries  that  are  parallel  to  the  cell's  rectilinear  lat¬ 
tice  unit  vectors  (e  and  y  in  Figure  3-15).  Such  boundaries  form  parallelograms. 
Radiation  from  irregularly  shaped  cells  is  treated  in  Section  3.4.3. 


Consider  a  coordinate  system  as  indicated  in  Figure  3-15.  With  respect 
to  this  system, 

r  ,  of  (3-163)  becomes  r'  where  (Section  3.2.4) 
oi  mn 


and 


-V  A  * 

r'  =  md  e  +  nd  y 

mn  e  y 


r  •  r'  =  md  s  +  nd  s 
mn  e  e  y  y 

and  the  field  point  e  and  y  directional  cosines  s  and  s  are 
r  e  y 


s  =  r  •  e 
e 


s  =  r-  y 

y 

The  of  (3-158)  become,  for  each  first  approximation  feedback  mode. 


(3-165) 

(3-166) 
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(3-168) 
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where 

a  =  a  e  +  a  v 
e  yJ 

Equation  (3-163)  then  becomes 
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(3-170) 
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(3-171) 


(3-172) 

(3-173) 


and  and  are  the  number  of  elements  along  constant  y  and  constant  e  lattice 
lines,  respectively,  in  the  cell.  These  numbers  are  =  6  and  =  7  for  the  Figure 
3-15  cell. 


cosines 


The  array  factor  can  be  expressed  in  terms  of  the  x  and  y  directional 


s  =  sin  0  cos  <t 
x 

s  =  sin  0  sin  <p 
y 


Since  (Section  3.2.4) 


e 
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Figure  3-16.  Definition  of  Three  Arrays,  Having  Closed  Form 
Array  Factors,  that  Combine  to  Form  Cell  Array 
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In  closed  form,  F  *  F(s  , 
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(3-177) 


(3-178) 


(3-179) 

(3-180) 


3.4.3  Irregular  Cells 

The  closed  form  expressions  for  the  cell  array  factor  developed  in  Sec¬ 
tion  3.4.2  assume  that  the  cell  boundaries  coincide  with  the  cell  lattice  vectors. 
This  restricts  cell  shapes  to  parallelograms.  In  general,  the  geometry  of  the  cell 
boundaries  may  not  coincide  with  the  rectilinear  structure.  An  example  is  shown  in 
Figure  3-16  (heavy  line).  A  closed  form  expression  for  the  array  factor  for  this  cell 
is  obtained  by  combining  several  array  factors.  Three  additional  arrays  are  indica¬ 
ted  in  Figure  3-16.  Array  1  is  "regular"  (bounded  by  a  parallelogram  with  sides  co¬ 
inciding  with  constant  lattice  coordinate  lines).  The  array  factor  for  Array  1, 
therefore,  is  given  by  (3-178).  Arrays  2  and  3  form  triangles  bounded  on  two  sides 
with  constant  lattice  coordinate  lines.  Along  each  third  side  the  y  coordinate  var¬ 
ies  linearly  with  e  coordinate.  A  closed  form  expression  for  this  array  type  is  de¬ 
rived  below.  The  cell  array  factor  then  is  given  by 


jkr.r.  jkr-r.  jkr-r, 

F  =  e  Fj^-e  F2+e  F3  (3-181) 

where  the  subscript  indicates  the  array  number.  The  phase  factors  kr*r^,  etc.,  are 
necessary  to  translate  the  corresponding  array  coordinate  origins,  references  for 
array  factors  F^ ,  etc.,  to  the  cell  coordinate  origin.  In  the  Figure  3-16  example. 


Array  factors  for  array  shapes  such  as  Array  2  and  Array  3  of  Figure 
the  form 
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n^O  - 

This  equation  is  similar  to  (3-171)  except  that  in  (3-182)  is  now 
m  (i.e.,  N  =  N^(m)).  A  linear  behavior,  expressed  by 


N  (m) 


e0  +  clm 


results  in  a  closed  form  expression  for  (3-182).  The  inner  summation 
becomes 
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This  results  in 
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which  reduces  further  to 
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3-16  are  of 

(3-182) 

a  function  of 

(3-183) 
in  (3-182) 

(3-184) 


(3-185) 


(3-186) 
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where 


<|/  =  (e^  +  <i»e>  /  2 


(3-187) 


and  iJj^  and  are  given  by  (3-179)  and  (3-180). 

Note  that  when  =  0,  Ny  is  a  constant  (=£q)  and  (3-186)  reduces  to  (3-178)  as 
expected. 


3.4.4 


Flat  Lens 


The  far-field  radiation  pattern  of  a  flat  lens  face  is  considered  here. 
Curved  lens  faces,  as  may  occur  from  structural  deployment  errors  or  from  severe 
temperature  gradients,  are  considered  in  Section  3.4.5. 

The  lens  patterns  are  with  respect  to  a  lens  coordinate  system,  termed 
"global,"  having  z  axis  directed  outward  from  the  lens  plane.  The  relationship  be¬ 
tween  the  global  coordinate  system  and  the  local  coordinate  system  for  the  v*1*1  cell 
is  shown  in  Figure  3-17.  The  g  subscript  refers  to  the  global  system  and  the  v  sub¬ 
script  to  the  cell  system.  The  vth  cell  is  positioned  rgv  from  the  global  system 
and  rotated  an  angle  <f>v  in  the  lens  plane  as  indicated. 

The  lens  pattern  is  obtained  by  combining  the  contributions  from  all  the 

cells.  Closed  form  cell  pattern  expressions  with  respect  to  local  cell  cool,  ■'nate 

systems  have  been  developed  in  Sections  3.4.2  and  3.4.3.  These  expressions  must  be 

modified  to  refer  to  the  global  coordinate  system.  The  array  factor  of  the  v^  cell 

F  (s  ,  s  )  is  given  by 
v  x  y 


F  (s  ,  s  ) 
v  x  y 


jkr-r 


8V  F(s  ,  s  ) 
x  y 


(3-188) 


where  F(s  ,  s  )  is  (3-178)  or  of  the  type  (3-181)  and 
x  y 


s  =  r  •  x 
x  v 


s  =  r-y 

y  J  v 


(3-189) 

(3-190) 


The  global  system  rectangular  components  of  the  field  point  unit  vector  r,  in  column 
vector  form,  are  given  by 


(3-191) 
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In  the  local  coordinate  system  of  the  v  cell,  r  has  rectangular  components  given  by 
r  •  x 

v 

■*•  =  f  U  v  /i  in'i\ 


!  R  r 
_  v  g 


(3-192) 


where  [R^]  is  the  coordinate  system  rotation  matrix.  For  flat  cells,  the  z  axes  of 
the  global  and  local  coordinate  systems  coincide  resulting  in 


cos  <k  sin  *  0 

v  v 

-sin  4>  cos  <(>  0 

V  V 


(3-193) 


Equation  (3-162)  for  the  element  factor  e  and  (3-178)  or  equations  of 
type  (3-181)  for  the  array  factor  F  employ  r  with  components  given  by  r^  in  deter¬ 
mining  the  vtl1  cell  contribution  to  the  lens  pattern. 

The  lens  fields  (3-164)  are  computed  from  the  6  and  ()>  components  of  the 
element  factor,  eg  and  e^ .  Here  6  and  4>  are  with  respect  to  the  global  system. 

Since  and  £  V  in  (3-162)  are  generally  known  in  local  system  coordinates,  a  deter¬ 
mination  of  en  and  e.  requires  that  C  and  $  be  expressed  in  local  coordinates.  In 
column  vector  form,  the  6  and  $  global  system  rectangular  components  are  respectively 


cos  6  cos  tj>  | 
cos  9  sin  <*>  j 
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(3-194) 
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(3-195) 


til 

The  corresponding  v  cell  rectangular  components  are 

9  =  [R  ]0 

v  vJ  g 

4»  =  [R  u 

v  v  g 

where,  for  flat  cells,  [R  ]  is  given  by  (3-193). 


(3-196) 

(3-197) 
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Figure  3-18.  Warped  Lens  Surface  S(x,y)  with  Quadrilateral  Approximation  Q(x,y) 


3.4.5  Nonflat  Lens 


The  pattern  expressions  —  (3-162),  (3-178),  and  (3-181)  —  assume  a  flat  lens 
array.  Departures  from  flatness  can  be  caused  by  the  adverse  space  envrionment  and 
possible  antenna  mechanical  malfunctions  or  imperfections.  An  example  of  the  former 
would  be  array  warping  caused  by  thermal  gradients  across  the  lens  surface  induced 
by  alternate  light  and  shadow  areas.  An  example  of  the  latter  would  be  an  imperfect 
mechanical  deployment  of  the  lens  in  space.  The  development  in  this  section  gener¬ 
alizes  work  done  by  Coffey  and  Carver.15 

The  total  lens  warping  is  represented  in  a  piecewise  manner  using  the  cells 
defined  in  Sections  3.4.2  and  3.4.3.  A  warped  array  is  shown  in  Figure  3-18.  If 
the  warp  is  smooth  and  slowly  varying,  as  in  thermal  expansion,  the  array  surface 
can  be  modeled  by  a  continuous  function,  .  Discontinuous  warps,  such  as  im¬ 

perfect  antenna  panel  unfolding,  can  be  modeled  by  piecewise  continuous  functions. 

If  the  array  edges  are  connected  by  straight  lines  (the  dotted  lines  in  Fi¬ 
gure  3-18),  the  resulting  surface  is  a  quadrilateral,  Q  (x^.y^),  that  approximates 
S(xv,yv)  linearly  in  the  arguments  x v  and  y  .  The  functional  form  for  Q(xv,yv)  is 

Q(x  ,y  )  =  b  +  b.x  +  b„y  +  b.x  y  (3-198) 

v  ■'v  0  1  v  2yv  3  vv 

here  the  coefficients  b„,  b, ,  b„ ,  and  b0  are  unknowns.  The  x  ,y  cross-term,  b,,  is 

0  1  2  3  v’-'v  ’  3 

a  measure  of  the  twist  in  the  warped  surface.  The  unknown  coefficients  of  Q(xv>yv) 

could  be  determined  by  choosing  four  "representative"  array  elements  locations 

(e.g.,  the  elements  closest  to  the  corners  of  the  array)  and  obtaining  four  linear 

equations  in  the  four  coefficient  unknowns.  Different  representative  elements  will 

lead  to  different  choices  of  Q(x  ,y  ). 

v  v’-'v 

There  should  be  little  variation  in  the  various  '  idrilaterals  for  small 
warps.  A  least-squares  analysis  could  be  performed  on  all  array  element  locations 
to  ootain  a  "best"  quadrilateral  approximation.  It  is  believed  that  such  an  analy¬ 
sis  will  not  be  necessary. 

A  simple  expression  for  the  array  factor  for  a  quadrilateral  surface  is  not 
presently  available.  Consequently,  a  further  approximation  must  be  made.  A  planar 
surface,  P(xv,yv),  of  the  form 

P(x  ,y  )  =  a_  +  a.x  +  a,y  (3-199) 

v  ■'v  0  1  v  2  v 

is  made  to  "best  fit"  the  quadrilateral  Q(x^,yv)  by  minimizing  the  integral 
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|q  -  P|dxdy 


(3-200) 


The  minimization  process  leads  to  the  following  optimum  choices  for  the  coefficients 
in  (3-199): 


a0  ~  b0 
al  =  bl 
a2  =  b2 


(3-201) 


The  plane,  P(x  ,y  ) ,  containing  the  optimum  coefficients  will,  in  general,  be  tilted 

with  respect  to  the  x^.y^  plane  (Figure  3-19).  The  tilt  is  measured  by  the  polar 

angles  6  and  6  .  The  equations  that  determine  the  tilt  angles  can  be  deduced  from 
x  y 

knowledge  of  P(xv,y^). 


tan  3  =  a, 

x  1 

tan  3  =  a. 

y  2 


(3-202) 


The  overall  result  is  that  the  warped  section  of  the  array  can  be  approxima¬ 
ted  by  a  lens  array  cell  tilted  in  the  local  unperturbed  cell  coordinate  system  an 
amount  meas-ired  by  the  tilt  angles.  Element  periodicity  is  assumed  to  still  exist 
on  the  tilted  cell  although  this  periodicity  may  differ  from  that  of  the  unperturbed 
cell.  Closed  form  array  factor  expressions  are  not  employed  if  the  elements  are  not 
periodically  located. 

The  pattern  functions  must  be  generalized  to  include  the  case  of  a  "tilted" 
planar  array.  Let  rj^  extend  from  the  origin  of  the  vtb  unperturbed  cell  local  co¬ 
ordinate  system  to  the  mntb  ground  screen  lattice  point  on  the  corresponding  tilted 
cell  (Figure  3-20).  Then 


-  ,  -*t 

r  =  a„z  +  r 
mn  0  v  mn 


(3-203) 


where 


rt=xx+yy+  (a.x  +  a_y  )z 
mn  mn  v  mn  v  1  mn  2  mn  v 


(3-204) 


and  x  ,  y  ,  and  z  are  the  cartesian  unit  vectors  of  the  v  unperturbed  cell  local 
v  v  v  ^t  th 

coordinate  system,  r^  locates  the  mn  lattice  point  in  the  tilted  cell  with  re¬ 
spect  to  a  local  tilted  cell  coordinate  origin. 


The  vector  r  defines  the  tilted  cell.  The  tilted  cell  lattice  vectors  and 
mn 

a  Suitable  tilted  cell  local  coordinate  system  can  be  determined  from  r  .  The  flat 

mn 
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Figure  3-20.  Unperturbed  Cell  and  Tilted  Cell  Approximation  to  a  Warp 


cell  radiation  pattern  expressions  developed  in  Sections  3.4.2  and  3.4.3  now  are  ap¬ 
plicable  to  the  tilted  cell.  The  contribution  of  the  tilted  cell  to  the  lens  pattern 
then  can  be  obtained  by  a  generalization  of  the  Section  3.4.4  development.  This  is 
carried  out  below. 

The  mnt'1  lattice  point  in  the  unperturbed  cell  is  defined  by 


r ' 
mn 


md  e  +  nd  y 
e  v  y  v 


(3-205) 


The  quantities  in  (3-205)  are  the  same  as  those  in  (3-165)  with  the  v  subscripts  im¬ 
plied  for  e  and  y  in  (3-165).  From  Section  3.2.4, 

r'  =  md  x  +  (nd  -  mAy)  y  (3-206) 

mn  x  v  y  3v 

The  x  and  y  in  (3-204)  are  the  coordinates  of  the  mnt^  lattice  point  in  the  un- 
mn  ■7mn 

perturbed  cell.  From  (3-206),  therefore, 

x  =  md  (3-207) 

mn  x 

y  =  nd  mAy  (3-208) 
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and  (3-204)  becomes 


r  =  md  x  +  (nd  -  mAy)y 
mn  xv  y  3  3v 


+  U-jmd^  +  a2(ndy  -  mAy)]  z v 


(3-209) 


The  y-directed  lattice  vector  of  the  tilted  cell  is  defined  by 


,t-t  -*t  -*t 

d  v  =  r  -  —  r 
y3  m,n+l  mn 


(3-210) 


where  d^  is  the  tilted  cell  y-directed  interelement  spacing  and  y1"  is  the  tilted  cell 
y-directed  lattice  unit  vector.  Similarly,  the  e  directed  lattice  vector  of  the  til¬ 
ted  cell  is  defined  by 

(3-211) 


,  t  - 1  -+t  -*t 

d  e  =  r  -  r 

e  m+l,n  mn 


with  similar  definitions  for  d1"  and  e*~ . 

e 

(3-211)  combine  to  yield 

d1"  =  d  /cos8 
y  y  y 


Equations  (3-202),  (3-209),  (3-210),  and 


^42/cos^  B  +  y2/cos2B  -  2dxAytanBxtan6 


cosB  y  +  sinB  z 
y  v  y  v 


(3-212) 

(3-213) 

(3-214) 
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e1"  =  —  [d  x  -  Ayy  +  (d  tan  8  -  Ay  tan  6  )z  ] 

^t  xv  JJv  x  x  }  y  v 

e 


(3-214) 


The  x  and  y  cartesian  unit  vectors  in  the  local  tilted  cell  coordinate  system  are 
obtained  readily  from 


xC  =  A/  /A- A 

-t  ~t  't 
z  =  x  x  y 


where 


t  «t  ,~t  -t.-t 
A  =  e  -  (y  -e  )y 


(3-216) 
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Equations  (3-212)  through  (3-218)  define  the  lattice  parameters  and  local  coordinate 
system  of  the  tilted  cell.  The  radiation  from  a  tilted  cell  then  can  be  determined 
from  the  expressions  of  Sections  3.4.2  and  3.4.3. 

A  tilted  cell  lattice  structure  may  differ  from  its  original,  unperturbed  lat¬ 
tice  structure;  e.g.,  an  unwarped  cell  with  a  rectangular  lattice  (Ay  =  0)  has,  after 
warping,  a  tilted  cell  approximation  with,  from  (3-213)  and  3-215), 


1 


Since  r  will  be  specified  with  respect  to  the  global  system  coordinates  (i.e.,  via 
the  column  vector  of  (3-191))  and,  since  r  with  respect  to  the  unperturbed  cell 
coordinates  is  r^  of  (3-192),  then,  with  respect  to  the  tilted  cell  coordinates,  the 
column  vector  representation  for  r  is 

rC  =  [RC]  [R]r  (3-224) 

v  g 

where  [R*1]  is  the  tilted  cell  rotation  matrix  given  by 


v 


•  x 

V 


-t 

X 


V 


•y 


V 


[RC] 


-t  - 

y  .x 

J  XT  11 


y  -y 

V  'v 


- 

y  • z 

J  XT  V 


(3-225) 
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The  expression  for  the  element  factor  (e)  for  a  tilted  cell  is  given  by  (3-162) 


with  r  having  components  given  by  r 


-t 


tilted  cell  coordinate  system,  a  determination  of  e0  and  e. 


Since  and  Z"  are  known  with  respect  to  the 


requires  that  0  and  4>  be 
expressed  in  tilted  cell  coordinates.  In  column  vector  form,  these  coordinates 
become 


=  [Rt] 


[R]  e_ 


(3-226) 


=  [R1]  [R]5  (3-227) 

where  0  and  $  are  defined  by  (3-194)  and  (3-195).  Equations  (3-196)  and  (3-197) 

8  8 

were  employed  in  deriving  (3-226)  and  (3-227). 


i 


SECTION  4 
IMPLEMENTATION 


The  SBR  RF  simulation  successive  approximation  method  software  is  subdivided 
into  two  sections,  based  upon  level  of  approximation.  One  section  implements  the 
First  Approximation  Method  (FAM) ,  which  makes  extensive  use  of  infinite  array  theory 
(Section  3.2).  The  second  implements  the  Higher  Approximation  Method  (HAM),  which 
employs  previous  approximation  results  in  a  conventional  isolated  element  moment 
method  approach  to  calculating  the  edge  effects  and  effects  of  nonperiodicities  on 
current  distributions  and  antenna  patterns  (Section  3.3). 

4.1  PROGRAM  FUNCTIONAL  DESCRIPTION 

4.1.1  First  Approximation  Method 

The  first  approximation  analysis  is  performed  for  each  cell.  Data  for  array 
a  are  read  in  and  the  element  data  are  calculated  in  subroutine  GEOMET  (Figure  4-1). 
GEOMET  is  called  again  to  calculate  array  b  element  data  after  array  b  data  are  read  in. 

The  program  loops  over  the  incident  plane  waves  from  1  to  NPWAVE.  An  inner 
loop  then  begins  over  the  feedback  modes.  The  feedback  mode  counters  are  calculated 
for  both  array  a  and  array  b.  The  phase  functions  for  the  total  impressed  field  are 
calculated  at  each  array  reference  element. 

Subroutine  ZMATX  is  called  with  array  b  data  to  calculate  the  array  b  impe¬ 
dance  matrix,  [Z],  for  an  infinite  array  above  an  infinite  perfect  ground  screen 
and  [ZF]  for  an  infinite  array  in  free  space.  Subroutine  SCAN  is  called  to  calcu¬ 
late  the  array  b  admittance  matrix  by  inverting  [Z]  and  calculating  the  active  impe¬ 
dance.  The  array  b  reference  element  current  for  unit  port  voltage  excitation  is 
saved.  The  reference  element  port  current  is  taken  from  the  LSB1"^  column  of  [Z], 
where  LSB  is  the  number  of  the  current  segment  containing  the  antenna  element  feed- 
point  for  array  b  elements. 

ZMATX  is  called  again  with  array  a  data  to  calculate  the  array  a  above  a 
ground  screen  impedance  matrix  [Z],  The  impedance  matrix  for  array  a  in  free  space 
[ZF]  is  also  calculated  in  ZMATX  and  is  written  to  a  file.  Matrix  [Z]  for  array  a 
is  sent  to  SCAN  to  calculate  the  admittance  matrix  for  array  a  and  the  active  impe¬ 
dance  for  array  a.  The  LSA1*1  column  of  [Z]  is  saved.  It  is  the  element  current  for 
unit  port  excitation. 

For  a  zero  feedback  loop,  matrix  [Z]  for  array  a  is  sent  to  SCIA  to  calcu¬ 
late  the  array  a  generalized  voltage  and  current  vectors.  The  reference  element 
port  current,  SCIAEX,  is  also  calculated. 


Array  a  and  array  b  reference  element  port  currents  are  then  calculated  for 
the  zero  feedback  loop  using  SCIAEX. 

If  there  is  feedback  and  if  the  number  of  feedback  modes  is  zero,  the  array  a 
impedance  matrix  for  free  space  [ZF]  is  read  from  the  file.  It  is  sent  to  SCAN  to 
calculate  the  array  a  reference  element  free  space  admittance  matrix  and  element 
active  impedance. 

The  array  a  impedance  matrix  for  free  space  [ZF]  and  the  array  a  active  impe¬ 
dance  are  sent  to  subroutine  MUTUAL  to  calculate  the  mutual  impedance  between  array 
b  and  array  a  through  the  imperfect  ground  plane.  MUTUAL  also  calculates  the  array 

a  generalized  voltage  and  current  vectors  arising  from  current  on  array  b.  The 

mutual  impedance  is  used  to  calculate  the  feedback  mode  array  a  and  array  b  refer¬ 
ences  element  port  currents. 

When  there  is  no  feedback,  the  call  to  SCAN  for  the  free  space  active  impe¬ 
dance  and  the  call  to  MUTUAL  are  bypassed. 

The  program  continues  by  calculating  the  reference  element  port  feedback  vol¬ 
tages.  For  feedback,  the  voltage  array  is  corrected  for  array  b  port  current  exci¬ 
tation.  For  no  feedback,  the  voltage  and  current  arrays  are  set  to  zero. 

Array  a  and  array  b  generalized  current  vectors  then  are  adjusted  for  the 
correct  array  port  voltages.  The  current  and  voltage  arrays  for  array  a  and  array 

b  are  written  to  file.  This  ends  the  feedback  loop. 

Subroutine  CELPAT  is  called  to  calculate  the  array  b  current  moments  for  the 
reference  elements  and  the  pattern  of  the  antenna  cell  under  analysis.  The  indivi¬ 
dual  cell  patterns  are  stored  on  a  file  for  integration  into  the  total  lens  pattern 

Array  a  reference  element  voltage  and  current  induced  for  each  incident  plane 
wave  are  written  to  file.  This  ends  the  plane  wave  loop. 
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Figure  4-1.  FAM  Flow  Chart 
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Figure  4-1.  FAM  Flow  Chart  (continued) 


Figure  4-1.  FAI1  Flow  Chart  (continued) 


4.1.2  Higher  Approximation  Method 

The  higher  approximation  method  (HAM)  calculates  corrections  to  the  first 
approximation  of  the  lens  patterns  and  current  amplitudes  to  account  for  edge  effects 
and  missing  elements,  or  any  nonperiodicities  in  the  array  lattice.  HAM  will  be  the 
implementation  of  the  theory  discussed  in  Section  3.3. 

HAM  begins  the  calculation  of  the  second  approximation  with  a  loop  over  the 
nonperiodicities  in  the  lens  (Figure  4-2).  For  each  nonperiodicity,  HAM  loops  over 
the  elements  in  the  neighborhood  of  the  nonperiodicity. 

The  physical  characteristics  of  the  radiating  elements  on  array  a  and  array 
b  are  read  in.  Next,  the  Z-parameters  of  the  module  are  input.  The  program  then 
calculates  the  generalized  admittance  row  matrices  for  array  a  and  array  b  elements 
above  a  perfect  ground  plane.  The  input  impedances  are  calculated  from  the  array 
a  and  array  b  elements'  admittances  which  are  extracted  from  the  admittance  row 
matrices. 

The  module  port  admittance  matrix  relating  voltages  induced  along  array  a 
and  array  b  elements  to  module  port  currents  is  calculated  using  the  impedance  ma¬ 
trices  and  the  Z-parameters. 

The  generalized  voltage  for  the  array  a  element  induced  by  the  first  approxi¬ 
mation  currents  on  array  b  is  calculated  using  the  first  approximation  prediction  of 
mutual  coupling  through  the  ground  screen. 

Next,  the  generalized  voltage  induced  on  array  a  element  by  excitation  field 
is  calculated. 

The  program  then  loops  over  the  neighboring  elements.  In  this  loop,  the 
generalized  voltage,  induced  on  the  array  a  element  by  first  approximation  currents 
on  the  neighboring  array  a  element  is  calculated.  This  generalized  voltage  is  sum¬ 
med  to  obtain  the  generalized  voltage.  The  generalized  voltage,  induced  on  array 
b  element  by  first  approximation  currents  on  the  neighboring  element  is  calculated 
and  summed  to  obtain  the  generalized  voltage.  This  ends  the  loop  over  neighboring 
elements . 

The  generalized  voltages  at  the  array  a  element  are  summed  to  obtain  the  vol¬ 
tage.  This  voltage  is  used  to  calculate  the  second  approximation  port  currents. 

The  module  port  currents  for  the  second  approximation  are  stored  and  the 
loop  over  the  elements  in  the  neighborhood  of  the  nonperiodicity  ends,  completing 
the  second  approximation  calculation. 


For  higher  approximations,  the  code  assumes  the  current  distributions  will 
change  by  a  complex  constant.  The  complex  constant  is  the  ratio  of  the  2nd  to  the 
1st  approximations  of  the  module's  port  currents.  Therefore,  the  higher  approxima¬ 
tion  begins  by  adjusting  the  voltage  induced  at  the  array  b  element  by  the  ratio  of 
the  2nd  to  the  1st  approximation  module  port  currents  on  the  array  a  side  of  the  lens. 

The  higher  approximation  generalized  voltages  are  calculated  by  looping  over 
all  the  elements  in  the  neighborhood  of  the  chosen  element.  The  second  approxima¬ 
tion  port  currents  for  the  module  are  used  to  adjust  the  1st  approximation  general¬ 
ized  current  vectors.  The  adjusted  currents  are  then  used  to  calculate  the  second 
approximation  generalized  voltage  vectors  for  the  chosen  element.  The  third  approx¬ 
imation  port  currents  are  then  calculated  using  the  second  approximation  general¬ 
ized  voltages.  The  third  approximation  port  currents  for  the  elements  are  stored 
and  used  to  calculate  new  current  moments.  The  current  moments  are  used  to  calcu¬ 
late  a  correction  field  pattern.  When  the  loop  over  the  elements  near  nonperiodi¬ 
cities  is  completed,  the  third  approximation  correction  field  pattern  is  also 
completed. 

The  higher  approximation  correction  calculations  are  complete  upon  completion 
of  the  loop  over  all  the  nonperiodicities  on  the  SBR  RF  array  lens. 
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Figure  4-2.  HAM  Flow  Chart 
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Figure  4-2.  HAM  Flow  Chart  (continued) 
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Figure  4-2.  HAM  Flow  Chart  (continued) 
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Figure  4-2.  HAM  Flow  Chart  (continued) 
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Figure  4-2.  HAM  Flow  Chart  (continued) 
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DATA  FLOW 


4.2.1  Data  Flow  for  FAM 

FAM  -  First  Approximation  Method  Main  Program 


Inputs 

ALPHAX,  ALPHAY 

AMP 

CELLID 

D 

DA 

DELTA(NBF,  / 
DELTAA(NBFA)  \ 

DX,  DY 

DXA,  DYA 

FREQ,  K 

LBA(NBFA, 3) , 
LB(NBF,3) 

L,  LSA,  LSB 

NA,  NB 

NBFA,  NBF 

NBFM 


Module  phase  shifts 

Amplitude  of  external  field  over  cell 

Alphanumeric  label  which  characterizes  the  antenna 
cell  under  analysis 

Height  of  array  b  above  the  ground  plane 
Height  of  array  a  above  the  ground  plane 

Arrays  of  lengths  of  wire  current  segments  {A} 

Array  b  interelement  spacing 

Array  a  interelement  spacing 

Radiation  frequency  and  wave  number 

Arrays  of  vectors  parallel  to  array  a  and  array  b 
wire  current  segments 

Number  of  the  current  segment  containing  antenna  element 
feedpoint  for  array  a  and  array  b  elements 

Number  of  elements  on  array  a  and  array  b  of  the  lens 

Number  of  basis  functions  which  describe  the  currents 
on  the  array  a  and  array  b  elements 

Number  of  feedback  modes  used  in  mutual  coupling  cal¬ 
culation 


NPWAVE 
PMAX,  QMAX 

PMAXA,  QMAXA 


Number  of  external  plane  waves  exciting  array  a 

Control  parameters  for  the  calculation  of  the  impedance 
matrix  for  array  b 

Control  parameters  for  the  calculation  of  the  impedance 
matrix  for  array  a 


PMAXM,  QMAXM 


Control  paremeters  for  the  calculation  of  the  array  b  to 
array  a  mutual  impedance 


(POLX,  POLY,  POLZ)  Polarization  vector  of  external  field 


RBCA(NBFA, 3) , 
RBC(NBF, 3) 


Array  of  position  vectors  R  pointing  to  the  center  of 
array  a  or  array  b  element  current  segments  from  array 
element  feedpoint 


(SX,  SY,  SZ) 


Direction  cosines  of  propagation  direction  of  external 
field 


TPERP ,  TPARA 


Perpendicular  and  parallel  ground  plane  transmission 
coefficients 


YDELA,  YDEL 


Array  a  and  array  b  skew  parameters 


ZMODAA,  ZMODAB  ( 
ZMODBA,  ZMODBB  \ 


Module  2-port  impedance  parameters 


Outputs 
IAAA(NBFA)  i 

IAAB(NBFA)  >  Array  a  reference  element  generalized  current  vectors 

IAEX(NBFA  ) 

IBBB(NBF)  Array  b  generalized  current  vector 

VAAA  j 

VAAB(NBF  ^  Array  a  voltages 

VAEX(NBF)  ' 

VBBB  Array  b  port  voltage 


GEOMET  -  Calculates  Array  Element  Data 
Inputs 


FREQ 

JCODE 

NFP 

NS 

NW 

R 

XB,  YB,  ZB,  j 

XE,  YE,  ZE  \ 

XF,  YF,  ZF 

Outputs 

DELTA (NBF) 

K 

LB (NBF, 3) 

NFS(5) 

NWS (NBF) 
RAD(NBF) 
RBC(NBF,3) 


Radiation  frequency 
Connector  code 
Number  of  feedpoint 
Number  of  segments 
Number  of  wires 
Input  radius 

Beginning  and  end  points 
Feedpoint  location 


Segment  length 
Wave  number 
Segment  unit  vector 

Array  containing  basis  function  number  for  each 
feed  segment 

Segment  wire  number 

Radius 

Vector  to  center  of  segment 


MUTUAL  -  Caluculates  the  Array  b  to  Array  a  Mutual  Coupling  Through  an  In 
finite  Imperfect  Ground  Plane 

Inputs 

D  Height  of  array  b  above  the  ground  plane 

DA  Height  of  array  a  above  the  ground  plane 


DELTA (NBF) 
DELTAA(NBFA) 
DX,  DY 
IBBB(NBF) 

K 

LB (NBF, 3) 
LBA(NBFA, 3) 


Array  b  lengths  of  wire  current  segments 
Array  a  l>-,.gths  of  wire  current  segments 
Array  b  interelement  spacing 

Array  of  current  segments  containing  array  b  feedpoint 
Wave  number 

Array  of  vectors  parallel  to  array  b  wire  current  segments 
Array  of  vectors  parallel  to  array  a  wire  current  segments 


Number  of  the  current  segment  containing  antenna  element 
feedpoint  for  array  b 

Number  of  basis  functions  for  array  b 

Number  of  basis  functions  for  array  a 

Control  parameters 

Flag  to  indicate  pulse  or  point  matching 

Array  of  position  vectors  pointing  to  center  of  array  b 
element  current  segments  from  array  b  element  feedpoint 

Array  of  position  vectors  pointing  to  center  of  array  a 
element  current  segments  from  array  a  element  feedpoint 

Array  b  phase  factors 

Skew  parameter  for  array  b 

Array  a  admittance  matrix 


Array  a  reference  element  generalized  current  vector 
due  to  VAAB  excitation 

Generalized  voltage  vector  on  array  a  due  to  penetration 
of  ground  screen  by  field  radiated  by  array  b  currents 

Mutual  impedance  between  array  b  and  array  a  reference 
elements  through  the  imperfect  ground  plane 

SCAN  -  Calculates  Array  Reference  Element  Admittance  Matrix  and  Active  Impe¬ 
dance 

Inputs 
L 

NBF 

Z(NBF,  NBF) 

Outputs 
Z(NBF,  NBF) 

ZSCAN 


The  number  of  the  current  segment  containing  antenna 
element  feedpoint  for  array  elements 

Number  of  basis  functions 

Impedance  matrix 

Inverted  impedance  matrix;  the  admittance  matrix 
Active  impedance 


LSB 

NBF 

NBFA 

PMAXM,  QMAXM 
PTPL 

RBC(NBF, 3) 

RBCA(NBFA, 3) 

TSXB,  TSYB 
YD  EL 

Z (NBFA,  NBFA) 
Outputs 

IAAB(NBF) 

VAAB (NBF) 

ZMUTAB 
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SCIA  -  Calculates  Array  Generalized  Voltage  and  Current  Vectors  and  Element 
Port  Current  for  Element  in  the  Incident  Field  ' 


Inputs 
DELTA (NBF) 
K 
L 

LB (NBF, 3) 


Array  of  lengths  {A}  of  wire  current  segments 
Wave  number 

Number  of  the  current  segment  containing  antenna  element 
feedpoint  for  array  element 

Array  of  vectors  [s? ]  parallel  to  array  wire  current  seg¬ 
ments 


NBF 

(POLX,  POLY,  POLZ) 

RBC(NBF, 3) 

Z(NBF,  NBF) 

OUTPUTS 

IAEX(NBF) 

SCIAEX 

VAEX(NBF) 


Number  of  basis  functions 
Polarization  vector  of  external  field 

Array  of  position  vectors  [R]  pointing  to  center  of  ar¬ 
ray  element  current  segments  from  array  element  feed- 
point 

Admittance  matrix 


Short  circuit  current  vector 

Short  circuit  source  current  for  reference  element 

Total  excitation  vector  (including  reflection)  for  re¬ 
ference  element 


ZMATX  -  Calculates  the  Impedance  Ilatrix 


Inputs 

D 

DELTA (NBF) 
DX,  DY 
K 

LB(NBF,3) 

NBF 
PTPL 
RAD  (NBF) 
RBC(NBF,3) 

WLEN 

YDEL 


Height  above  the  ground  plane 

Array  of  lengths  (A)  of  wire  current  segments 

Interelement  spacing 

Wave  number 

Array  of  unit  vectors  [it]  parallel  to  wire  current  seg¬ 
ments 

Number  of  basis  functions 

Flag  to  indicate  0  =  point  matching  1  =  pulse  matching 
Radius 

Array  of  position  vectors  [R]  pointing  to  center  of  ele¬ 
ment  current  segments. 

Wave  length 

Skew  parameter 
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Outputs 


Z  (NBF ,  NBF) 


IV (NBF ,  NBF) 


Impedance  matrix  for  an  array  above  a  perfect  ground 


Impedance  matrix  for  an  array  in  free  space 


CELLPAT  -  Cell  Pattern  Calculations 


Inputs 


RGV ( 3 ) 


NY(30) 

NX(30) 

N0(30) 

IAS (30) 
RREFC(30,3) 
DX (30) 

DY (30) 
DELTAY(30) 
RBC(NBF,3) 
IBBB(NBF) 

Outputs 

THETA 

PHI 

ETHETA 


4.2.2  Data  Flow  for  HAM 


Main  Program 


Inputs 


ALPHAX,  ALPHAY 


CELLID 


Global  coordinates  of  cell 

Rotation  of  cell  relative  to  GCS 

Rotation  of  element  relative  to  CCS 

Number  of  arrays  within  cell 

Number  of  rows  in  y  direction  for  each  array 

Number  of  columns  in  x  direction  for  each  array 

Number  of  elements  in  first  column  of  irregular  array 

Add-subtract  factor  for  each  array 

Coordinates  of  reference  element  of  each  array 

x  separation  for  each  array 

y  separation  for  each  array 

Skew  in  y  direction  for  each  array 

Array  position  vectors  for  segments  of  array  a  element 
Generalized  current  vector  for  array  b  elements 


Polar  angles  to  field  points 
Azimuth  to  field  points 
Theta  component  of  E  field 
Phi  component  of  E  field 


Module  phase  shifts 

Amplitude  of  external  field  over  cell 

Alphanumeric  label  which  characterizes  the  antenna  cell 
under  analysis 


D 

DA 


DX,  DY 
DXA,  DYA 


FREQ,  K 


IAAA(NBFA) 

IAAB(NBFA) 

IAEX(NBFA) 


1 


s 


IBBB(NBF) 


LBA(NBFA,3), 
LB (NBF, 3) 


L,  LSA,  LSB 
NA,  NB 
NBFA,  NBF 


NFBM 


Height  of  array  b  above  the  ground  plane 
Height  of  array  a  above  the  ground  plane 
Array  b  interelement  spacing 
Array  a  interelement  spacing 
Radiation  frequency  and  wave  number 

Array  a  reference  element  generalized  current  vectors 
for  first  approximation 

Array  b  generalized  current  vector  for  first  approxi¬ 
mation 

—b 

Arrays  of  vectors  l  parallel  to  array  a  and  array  b 
wire  segments 

Number  of  the  current  segment  containing  antenna  element 
feedpoint  for  array  a  and  array  b  elements 

Number  of  elements  on  array  a  and  array  b  of  the  lens 

Number  of  basis  functions  which  describe  the  currents  on 
the  array  a  and  array  b  elements 

Number  of  feedback  modes  used  in  mutual  coupling  calcu¬ 
lation 


(POLX,  POLY,  POLZ)  Polarization  vector  of  external  field 


RBCA(NBFA, 3) , 
RBC(NBF, 3) 


Array  of  position  vectors  R  pointing  to  the  center  of 
array  a  or  array  b  element  current  segments  from  array 
element  feedpoint 


(SX,  SY,  SZ) 


Direction  cosines  of  propagation  direction  of  external 
field 


VAAA  i 

VAAB(NBF)  [ 
VAEX(NBF)  ) 

VBBB 

(XLOC,  YLOC,  ZLOC) 


Array  a  voltages 


Array  b  port  voltage 

Location  of  mn^  element  on  cell 


YDELA,  YDEL 


Array  a  and  array  b  skew  parameters 


ZMODAA ,  ZMODAB,  / 
ZMODBA,  ZMODBB  i 


Module  2-port  impedance  parameters 


Outputs  (Subscript j  denotes  order  of  approximation) 

IAAAj  (NBFA) 

IBBBj  (NBFA) 


AEOj 

AE4>j 

VAAAj 

VBBBj 


Array  a  generalized  current  vector 
Array  b  generalized  current  vector 


Electric  field  corrections  provided  by  higher  order 
approximations  of  voltages  and  currents 

Array  a  port  voltage  for  00  element 

Array  b  port  voltage  for  00  element 
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SECTION  5 


RESULTS  -  FIRST  APPROXIMATION  METHOD  CURRENTS  AND  ACTIVE  IMPEDANCE 

Simulation  results  obtained  to  date  have  been  limited  to  infinite  array 
analysis.  Active  impedance  and  current  distributions  have  been  compared  with  pub¬ 
lished  results  on  planar,  rectangular  lattice,  and  dipole  arrays . 1 1 ’ 1 2 > 1 3  Simula¬ 
tion  results  are  presently  being  compared  with  moment  method  analyses  of  large,  tri¬ 
angular,  spaced-dipole  arrays.  The  center  element  impedance  and  current  have  been 
compared  with  the  first  approximation  method  (FAM) .  Several  of  these  comparisons 
are  presented  in  this  section.  The  effect  of  varying  the  orientation  of  the  dipoles 
in  an  infinite  triangular  array  is  presented,  as  is  a  design  of  a  single-port  turn¬ 
stile  radiator  for  use  in  an  infinite  triangular  array.  These  results  demonstrate 
the  versatility  of  the  simulator.  Convergence  curves  are  given  in  terms  of  minimum 
number  of  required  plane  wave  modes  (truncated  doubly  infinite  summation  in  (3-16)) 
and  number  of  subsections  (order  of  [ZSS]). 

5.1  INFINITE  ARRAY  COMPARISONS 

Stark11  analyzed  an  infinite  rectangular  array  of  thin  A/2  strip  dipoles  over  a 
perfect  ground  screen.  Active  impedances  were  computed.  A  single  sinusoidal  current 
distribution  on  each  element  was  assumed  since  this  is  exact  for  vanishingly  thin 
wires  and  since  such  an  approximation  proved  valid  in  calculating  impedances  of  iso¬ 
lated  dipoles  not  exceeding  A/2  in  length.  These  results  compared  well  with  Carter's 
large  finite  array  analysis14  for  active  impedances  near  the  center  of  the  array; 
however.  Carter  also  assumed  the  same  single  sinusoidal  current  distribution  per 
element . 

VanKoughnett  and  Yen12  observed  that  since  the  dipoles  were  almost  touching 
in  the  collinear  direction  the  current  distribution  may  significantly  deviate  from 
a  single  sinusoid  even  for  very  thin  dipoles.  They  analyzed  the  same  array  in  a  man¬ 
ner  which  permitted  computing  the  current  distribution  to  any  desired  degree  of  ac¬ 
curacy.  The  resulting  active  impedances  displayed  a  variation  with  scan  that  was 
similar  to  Stark's  results;  however,  the  variation  was  considerably  more  pronounced 
and  the  broadside  impedance  values  significantly  differed. 

VanKoughnett  and  Yen's  method  is  applicable  only  to  elements  with  small  spacings 
in  the  collinear  direction.  Chang13  removed  this  restriction  by  applying  a  "five- 
term  expansion"  as  well  as  a  multiple  mode  sinusoidal  expansion  to  the  rectangular 
infinite  array  of  dipoles. 

The  results  from  the  above  papers  selected  for  comparison  with  the  first  approx¬ 
imation  method  (FAM)  were  obtained  with  the  configuration  of  Figure  5-1  where 
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a  =  dipole  radius 
d  =  dipole  length 

h  =  array  height  above  ground  screen 
d^  =  periodicity  in  x  direction 

d^  =  periodicity  in  y  direction 

The  excitations  were  unit  amplitude  voltages  progressively  phased  for  radiation  in  a 
specified  direction.  The  dimensions  for  three  problems  that  were  considered  for  '•his 
configuration  are  indicated  in  Table  5—1.  The  voltage  sources  were  phased  for  a  vari¬ 
ety  of  E-plane,  H-plane,  and  diagonal  plane  scan  angles.  Both  current  distributions 
and  active  impedances  were  computed.  The  results  from  Problems  1  and  2  apply  to 
Chang's  results.  The  results  from  Problem  3  apply  to  VanKoughnett  and  Yen's  results. 

The  zero  scan  (broadside  radiation)  reference  element  current  distribution  for 

Problem  1  is  shown  in  Figure  5-2.  Six  FAM  results  appear.  Each  result  applies  to 

one  of  six  combinations  of  maximum  plane  wave  expansion  numbers  p  and  q  and 

max  hnax 

number  of  segments  Ns  used  in  computing  the  infinite  array  generalized  impedance  ma¬ 
trix  [Zgs]  of  Section  3.2.2.  Equation  (3-16)  shows  the  doubly  infinite  summation 
in  p  and  q  needed  to  compute  the  ijth  element  of  the  N  x  N  matrix  [Zgs].  These 


summations  are  truncated  according  to  Ip  lip  and  Iql-Sq  .  The  p 

1  1  ‘max  1  ’ 1  hnax  *n 

Ng  for  the  FAM  results  in  Figure  5-2  are  listed  below. 


q  ,  and 
^max 


Curve  No. 

N 

p 

q 

s 

max 

’max 

1 

11 

15 

40 

2 

11 

15 

30 

3 

11 

15 

20 

4 

21 

15 

20 

5 

21 

15 

40 

6 

21 

15 

30 

required  p 

max 

and  q 

max 

increases  with 

increasing  d^, 

Chang's  results  (Problems  1  and  2)  d  is  almost  5  times  d  .  Therefore,  p  was 

x  y  rmax 

chosen  smaller  than  qmax-  The  Problem  1  current  distribution  results  (Figure  5-2) 
appear  converged  for  21  segments  per  half  wavelength.  The  curves  for  11  segments 
have  not  yet  converged.  Curves  1  and  2  are  further  from  Chang's  results  than  curve 
3  even  though  the  value  of  q  is  higher  for  curves  1  and  2  than  for  3.  Additional 
convergence  results  are  given  in  Section  5.5. 


5-2 


The  principal  plane  scan  active  admittances  for  Problem  1  are  plotted  in  Figure 

5-3.  For  these  results  p  =  20,  q  =70  and  N  =  21. 

rmax  ’max  s 

The  element  current  for  the  full  wavelength  dipole  array  of  Problem  2,  scanned 

for  broadside  radiation  is  shown  in  Figure  5-4.  The  FAM  results  used  p  =  15, 

rmax  ’ 

=  50  and  Ng  =  21.  As  with  Problem  1,  close  agreement  between  FAM  and  Chang's 
method  is  indicated. 

The  E-plane  and  H-plane  active  impedances  for  Problem  3  are  compared  with 
VanKoughnett  and  Yen  in  Figures  5-5  (resistance)  and  5-6  (reactance).  The  resistance 
is  in  close  agreement  throughout  the  0°  (broadside)  to  45°  scan  range.  The  agreement 
in  reactance  is  within  approximately  15%  for  near  broadside  scan.  This  agreement  de¬ 
teriorates  for  larger  scan  angles. 
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a.  Front  View 
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b.  Side  View 


Figure  5-1.  Infinite  Planar  Dipole  Array  Above  a  Perfect  Ground  Screen 
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TABLE  5-1 


INFINITE  ARRAY  PROBLEMS 


a 

d 

h 

d 

x 

d 

y 


Problem  1 

Problem  2 

Problem  3 

0.007022 

0.007022 

0.00796 

0.5 

1. 

0.5 

0.25 

0.25 

0.25 

0.25 

0.25 

0.5 

1.2 

2.2 

0.5 

All  dimensions  in  wavelengths 


Figure  5-2.  Element  Current  for  Problem  1  of  Table  5-1 


Figure  5-4.  Element  Current  for  Problem  2  of  Table  5-1 
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Figure  5-5.  Active  Resistance  for  Problem  3  of  Table 


Figure  5-6.  Active  Reactance  for  Problem  3  of  Table  5-1 


5.2  FINITE  ARRAY  COMPARISONS 


Thin  wire  moment  method  computer  codes  are  applicable  to  finite  arrays  of  arbi¬ 
trarily  bent  wire  radiators  over  infinite  perfectly  conducting  ground  screens.  Two 
of  these  codes  GEMACS18  and  AMP19  are  available  as  an  aid  in  validating  the  simulator. 

The  results  in  Section  5.1  apply  only  to  rectangular  lattices.  Since  triangular 
lattices  are  of  particular  interest,  the  configuration  chosen  for  comparing  FAM  with 
finite  array  moment  methods  is  that  of  Figure  5-7.  A  19-element  dipole  array  is  in¬ 
dicated.  Each  dipole  is  centrally  excited  with  progressively  phased  voltage  sources. 
The  parameter  values  are 

d  =  0.51 

a  =  0.011 

d  =  0.7151 
x 

d  =  0.8261 

y 

A  =  0.4131 

y 

h  =  0.251 


Consider  the  standard  spherical  coordinates  r,  6,  if  with  respect  to  the  right  hand 

x,  y,  z  coordinates  indicated  in  Figure  5-7  where  z  is  outward  normal  to  the  array 

plane  and  the  origin  is  at  the  reference  element  (m  =  n  =  o)  center.  The  phasing 

of  the  V  voltage  sources  can  be  related  to  main  beam  angular  directions  6  ,  if 
mn  °  b  a’  a 

(scan  angles)  by  letting  r^  denote  the  mnth  element  center  position  vector.  Then, 

with  reference  to  Figure  5-7  and  the  discussion  in  Section  3.2.4,  (3-95)  results  in 


-jkr*r  -ik(md  a  +  (nd  -  mAy)a  ) 

,,  mn  J  xx  y  y 

V  =  e  =  e  J  } 

mn 


(5-1) 


where  ct  and 

X 

a 

y 

are  the  scan  directional  cosines 

a  =  cos 

sin  6 

X 

a 

a 

a  =  sin 

$ 

sin  0 

y 

a 

a 

and  r  is  the  unit  vector  for  the  main  beam  field  point  position  vector.  Equation 
(5-1)  provides  the  element  excitations  for  any  scan  direction.  The  "principal" 
scans  are 

H-plane:  -  0  (<fa  =  0) 
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E-plane:  ax  ~  0  (<t  =  tt/ 2 ) 

Diagonal  plane:  ax  =  ay  (4>a  =  */4) 

FAM  treats  the  Figure  5-7  problem  as  an  infinite  array.  The  active  impedances 
for  the  principal  scan  planes  for  this  problem,  as  computed  by  FAM,  are  shown  in  Fig¬ 
ures  5-8  (resistance)  and  5-9  (reactance).  The  scan  angle  is  9  .  The  computational 

variables  were  N  =  23  and  p  =  q  =  60.  The  discontinuities  in  these  curves 
s  rmax  ’max 

occur  at  grating  lobe  singularities;  i.e.,  at  scan  angles  for  which  there  exist  in¬ 
tegers  p  and  q  such  that  q^ ,  given  by  (3-97),  vanishes.  (A  discussion  on  grating 
lobe  singularities  appears  in  the  paragraph  following  (3-58)  and  in  the  paragraph 
containing  (3-61).)  The  H-plane  grating  lobe  singularity  occurs  at  6^  =  23.49°,  the 
E-plane  singularity  occurs  at  6^  =  29.22°,  and  the  diagonal  plane  singularity  occurs 
at  =  59.75°.  The  locations  of  grating  lobe  singularities  depend  only  upon  lat¬ 
tice  spacings. 

The  broadside  scan  current  on  the  reference  element  (m  =  n  =  0)  of  the  Figure 
5-7  finite  array  as  computed  by  the  GEMACS  moment  method  code  is  shown  in  Figure  5-10. 
Also  shown  are  the  GEMACS  current  with  the  array  extended  to  include  37  dipoles 
(same  interelement  spacings)  and  the  infinite  array  FAM  current. 
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Figure  5-7.  Finite  triangular  lattice  dipole  array  a  height  h  above 
infinite  ground  screen.  Identifying  m,n  indices  are  in 
dicated.  A  voltage  source  centrally  excites  each 
dipole.  mn 


At  point*  where  •ore  than  one  curve  *eet>  letter*  occurring  earliest  alphabetical  lu  are  suppressed 
e*p.»  a  point  coaaon  ta  curves  A»C»  and  F  is  marked  F. 


Figure  5-8.  Active  Resistance  vs.  Scan  Angle  for  Figure  3-7  Problem  Treated  as  an  Infinite  Array 
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5.3  ELEMENT  ROTATION  IN  THE  ARRAY  PLANE 


The  effects  of  rotating  all  dipole  array  elements  in  the  array  plane  (about  the 
z  axis)  are  of  interest.  The  reference  element  currents  and  active  impedances  were 
obtained  using  FAM  for  rotation  angles  (Figure  5-11)  =  0°,  15°,  30°,  45°,  60°,  and 

90°.  This  was  done  at  both  broadside  scan  and  0  =  20°,  6  =  0°  scan.  Two  lattices 

a  a 

were  considered:  the  Figure  5-7  triangular  lattice  and  a  d^  =  =  0.8X  rectangular 

lattice.  The  active  impedances  are  given  in  Table  5-2.  The  reference  element  cur¬ 
rents  are  plotted  in  Figures  5-12  to  5-19.  These  results  indicate  greater  element 
rotation  induced  impedance  and  current  variations  for  the  20  degree  scan  than  for 
broadside  and  for  the  rectangular  lattice  than  for  triangular. 
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TABLE  5-2 


r 


ACTIVE  IMPEDANCES  FOP  ROTATED  ELEMENT  ARRAYS 


Lattice 

Scan  Angle  6^ (Degrees)  ($a=0°) 

Rotation  Angle  iKDegrees) 

Impedance (Ohms) 

Triangular 

0 

0 

78.87  +  j  8.614 

15 

78.4  +  j  7.773 

30 

77.99  +  j  7.147 

45 

78.41  +  j  7.803 

60 

78.87  +  j  8.607 

90 

77.86  +  j  6.806 

20 

0 

99.48  +  j41. 17 

15 

96.74  +  j39.13 

30 

90.9  +  J35.51 

45 

85.35  +  j  33. 52 

60 

81.56  +  j34.5 

90 

78.73  +  j37.95 

Rectangular 

0 

0 

78.71  +  j24.42 

15 

77.36  +  j21. 7 

30 

75.2  +  j 17. 19 

45 

74.26  +  jl5. 17 

60 

75.19  +  jl7. 18 

90 

78.72  +  j24.44 

20 

0 

166.6  +  j 12. 52 

15 

158.7  +  j 15. 66 

30 

140.1  +  j23. 63 

45 

118.3  +  j32.58 

60 

99.78  +  j39.53 

90 

83.98  +  j44.4 

o 


Figure  5-12.  Real  Part  of  Reference  Dipole  Current  for  Triangular  Lattice  (Figure  5-7) 
for  Several  Dipole  Orientations  and  Broadside  Scan  (0  =0°) 


point  common  to  curves  A»C»  and  F  is  marked 


At  points  uhere  more  than  one  curve  meet,  letters  occurring  earliest  alphabetically  are  suppressed 
8.9.1  a  point  common  to  curves  A»C»  and  F  is  marked  F* 


At  points  whir#  more  than  ort#  curve  meet,  letters  occurring  earliest  a  1  phabe  t  i  ca  1 1  y  are  suppressed 
e.S«»  a  point  common  to  curves  AtCr  and  F  is  marked  F. 


5.4  TURNSTILE  RADIATORS 

The  first  approximation  method  was  applied  to  the  design  of  a  turnstile  radiator 
(crossed  dipoles  excited  for  circular  polarization)  for  use  as  an  array  element.  Re¬ 
actance  loads  were  determined  to  effect  the  90  degree  phasing  between  horizontal  and 
vertical  arms.  This  design  avoids  having  to  excite  the  radiator  with  the  convention¬ 
al  two-port  source.  The  array  element  is  shown  in  Figure  5-20.  The  array  lattice 
was  that  of  Figure  5-7.  Only  broadside  array  phasing  was  considered,  resulting  in 
symmetric  current  about  the  axes  of  both  dipoles.  The  need  to  model  the  short  exci¬ 
ted  arm  was  avoided  by  passing  the  source  through  the  junctions.  This  resulted  in 
the  equivalent  problem  of  Figure  5-21. 

The  reference  element  current  as  computed  by  FAM,  with  the  complex  impedance 

loads  Z  =  Z  =  0  is  shown  in  Figures  5-22  and  5-23. 

H  V 

The  current  under  non-zero  load  conditions  is  easily  determined  by  FAM  as  dis¬ 
cussed  in  Section  3.2.7.  Also,  by  considering  the  radiator  as  a  multiport  -  one 
port  for  each  load  and  excitation  -  and  using  FAM  to  obtain  the  port  parameters  re¬ 
lating  these  values  of  Z^  and  Zy  were  determined  such  that  the  currents  at  the  load 
ports  were  equal  in  amplitude  and  in  phase  quadrature  between  horizontal  and  vertical 
arms.  The  resulting  complex  loads  are 

ZH  =  0.  +  j  32  (5-2) 

Z  =  -0.2  -  j48.8  (5-3) 

Since,  for  energy  considerations,  purely  reactive  loads  are  desirable  the  small  real 
part  of  Zy  was  set  to  zero  prior  to  computing  the  reference  element  current.  This 
current  is  plotted  in  Figures  5-24  and  5-25.  The  active  impedance  for  this  radia¬ 
ting  element  in  the  Figure  5-7  lattice  infinite  array  is  74.4  +  jl0.9  ohms. 
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Figure  5-20.  Loaded  Crossed  Dipoles  with  Single  Feed  Port 


igure  5-21.  Loaded  Crossed  Dipoles  with  Equivalent  Feed 
This  Model  is  Equivalent  to  the  Figure  5-20 
Single  Feed  Port  Radiator 
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where  *ore  than  one  curve  meet  letters  occurring  earliest  a  1 phahet i ca 1 1 *  are  suppressed 
int  c otfctuun  to  curves  A»C»  and  F  is  parked  F. 


POSITION  RELATIVE  TO  C ENTER  OK  CROSSED  HI  POl.ES<WAVE  LENGTHS) 

Figure  5-25.  Phase  of  Current  on  Crossed  Dipole  Array  Reference 
Element  with  Loading  for  Circular  Polarization 


5.5  CONVERGENCE 


Each  element  of  the  generalized  impedance  matrix  [Z8S],  given  by  (3-61),  in¬ 
volves  a  doubly  infinite  summation  in  integers  p  and  q.  Adequate  truncation  of  this 

summation  depends  on  its  rate  of  convergence.  As  discussed  in  Section  5.1,  p  and 

max 

q  must  be  determined  such  that  the  exclusion  of  all  p  and  q  for  which  Ip  I  >  p 
max  r  n  u  1  ‘max 

and  [ q |  >  q  in  the  summation  provides  an  accurate  result.  Each  pq  term  in  the 
summation  represents  a  plane  wave.  Beyond  the  first  few  values  of  |p|  +  |q|  these 
waves  are  evanescent;  they  exponentially  decay  away  from  the  array  plane.  The  sum¬ 
mation  term  corresponding  to  an  evanescent  wave  is  imaginary  since  the  energy  impart¬ 
ed  to  these  waves  is  entirely  reactive.  The  principal  parameters  which  impact  Pmax 
and  qmax  for  wire  radiators  lying  in  the  array  plane  are  d^/a  and  d^/a.  The  larger 

these  ratios  of  interelement  spacing  to  wire  radius  the  slower  the  convergence  of 
the  summation. 

The  reactance  of  a  "self- impedance"  (diagonal)  element  of  [Z8S]  was  computed 

for  values  of  p  =  q  ^  150  and  several  values  of  d  /a  =  d  /a  (square  lattice), 
max  ’max  x  y 

Convergence  rates  are  apparent  from  these  results  by  plotting  (Figures  5-26  to  5-29) 


e(M)  = 


,X(M)  -  X(150) 
1  X(150) 


(5-4) 


Og 

versus  M  where  X(M)  =  imaginary  part  of  a  diagonal  element  of  [Z&  ]  determined  with 
p  =  q  =  M.  Figures  5-26  to  5-29  each  correspond  to  a  different  wire  radius 
and  each  display  six  curves.  Each  curve  corresponds  to  a  different  interelement 
spacing.  Curve  C  in  Figure  5-26  and  Curve  A  in  Figure  5-28  correspond  to  the  same 
ratio  50,  of  interelement  spacing  to  radius.  As  expected,  these  two  curves  exhibit 
similar  convergence  rates. 

og  2 

The  size  of  [Z6  ]  is  Ng  where  Ng  is  the  number  of  segments  dividing  a  wire  ra¬ 
diator.  Convergence  of  the  reference  element  wire  current  with  increasing  Ng  is 
shown  in  Figures  5-30  to  5-33.  Figures  5-30  and  5-31  correspond  to  the  Figure  5-7 
triangular  lattice  dipole  array  and  p  =  q  =  60.  Figures  5-32  and  5-33  corre- 

TT13.X  ITlelX 

spond  to  the  Figure  5-7  array  with  the  exception  that  the  wire  radius  was  reduced  to 

a  =  0.005A  and  p  =  q  =  80.  Convergence  is  quite  slow  in  both  cases.  It  was 
max  max 

pointed  out20  that  the  use  of  pulse  expansion  functions  in  the  moment  method  formula¬ 
tion  in  conjunction  with  an  exact  representation  for  the  E-field  operator  is  likely 
the  cause  of  the  slow  convergence.  The  convergence  is  expected  to  be  significantly 


improved  by  computing  the  E-field  via  a  partial  difference  approximation, 
mately  seven  subsections  per  wavelength  is  expected  to  then  be  adequate, 
fication  is  presently  being  implemented. 


Approxi- 
This  modi- 


At  points  whpre  more  than  one  curve  meet,  letters  occurring  earliest  alphabetically  are  suppressed* 
e**t*»  a  point  common  to  curves  A»C*  arid  F  is  marked  F. 


Figure  5-26#  Convergence  of  Diagonal  Element  of  [Zg  ]  for  Wire  Radius  *  0.01A 


At  points  where  nore  than  one  curve  meet^  letters  occurring  earliest  alphabetically  are  suppressed* 
e.a.t  a  point  common  to  curves  A*C*  and  F  is  marked  F. 


Figure  5-27.  Convergence  of  Diagonal  Element  of  [Z  ]  for  Wire  Radius  =  0.007  A 


At  points  where  wore  than  one  curve  meet,  letters  occur  rina  earliest  a } phabe t i ca 1 1 v  are  suppressed 
e.*.»  a  point  cow won  to  curves  A*C*  and  F  is  marked  F, 


Figure  5-28.  Convergence  of  Diagonal  Element  of  [Zs  ]  for  Wire  Radius  =  0.005 A 


Figure  5-29.  Convergence  of  Diagonal  Element  of  [Ze  ]  for  Wire  Radius  -  0.001X 
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Figure  5-30.  Convergence  of  Figure  5-7  Reference  Element  current 
(Real  Part)  with  Respect  to  Number  of  Segments 
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APPENDIX  A 


FIELD  FROM  AN  INFINITE  ARRAY  OF  INFINITESIMAL  CURRENT  DIPOLES 


This  appendix  contains  the  derivation  of  the  equations  necessary  to  determine 
the  radiation  from  an  infinite,  rectilinear  array  of  infinitesimal  current  elements. 
All  elements  are  identical  in  orientation,  progressively  phased,  and  of  constant  cur¬ 
rent  amplitude  (Figure  A-l).  The  derivation  is  based  on  the  work  of  Munk,  Burrell, 
and  Komban.1’6 


The  vector  potential  for  an  infinite  array  can  be  written  as  a  sum  over  the 
vector  potentials  of  the  individual  elements 


dA  = 


E  E 


dA 


q=-o 


PS 


(A-l) 


where  the  ordered  integer  pair  (p,q)  uniquely  locates  an  element  in  the  array.  The 
individual  element  vector  potentials  have  the  form 


I 


pq 


(A-2) 


where  i  is  the  element  orientation,  p  is  the  permeability  of  free  space,  df.  is  the 

tVi 

infinitesimal  diode  length  ,  I  is  the  current  on  the  pq  element, R  is  the  dis- 

•  th  pq  pq 

tance  from  the  pq  element  to  the  field  point.  The  propagation  constant,  y,  is  as¬ 
sumed  complex  to  add  generality  to  the  argument. 


H 


The  distance  R 

pq 


is  given  by 


R 

pq 


(A-3) 


where  r  points  to  the  field  point  and  points  to  the  pq^  current  element.  These 
vectors  have  the  following  rectangular  coordinate  representations: 


r  =  xx  +  yy  +  zz 

-V.  *  « 

r  =  x  x  +  y  y 

pq  pq  pq 


(A-4) 

(A-5) 


where  the  array  is  assumed,  for  simplicity,  to  lie  in  the  xy  plane  (Figure  A-l). 

For  the  case  of  a  rectilinear  lattice  on  the  array,  the  coordinates  of  the  pq1*1  ele¬ 
ment  satisfy 

x  =  pd  (A-6) 

pq  * 

y  ”  qd  +  pAy  (A-7) 

pq  y 

A-l 


where  and  d  are  interelement  separations  and  Ay  is  the  array  skew  paremeter. 

(See  Figure  3-7  for  an  illustration  of  these  parameters.)  The  skew  parameter  will 
be  zero  for  rectangular  lattices.  The  distance  R  now  can  be  written  as 

pq 

R2  =  (pd  -  x)2  +  (qd  +  pAy  -  y)2  +  z2  (A-8) 

pq  x  ^  y  r  J 

Due  to  the  constant  amplitude  and  progressive  phasing  of  the  element  currents,  the 

current  I  of  the  pqt'1  element  can  be  written  in  terms  of  a  reference  current  as 
pq  00 

-j kr '  -s 

I  =  I  e  Pq  (A-9) 

pq  00 

where  k  is  the  wave  number  (k  =  imaginary  part  of  y)  and  s  is  the  phase  unit  vector. 


s=sx  +  sy  +  sz 
x  y  z 


The  current  X  in  (A-9)  now  assumes  the  form 

pq 


(A-10) 


-3kpVx  ~  3k(9dy  +  PAy>sy 


1  =x  e  1  J 

pq  00 

The  array  vector  potential  can  be  written  as 


(A-ll) 


dA  =  l 


PlQ0<^g'  V 

4tt 


-jkpds  -  jkpAys 
e  x  x  y  dA 


V  E 


-lkqd  s  -yR 

y  r  P 

e  e 


(A-12) 


(A-13) 


and  represents  the  vector  potential  for  a  row  of  elements.  The  distance  R  is 

pq 

written  for  convenience  as 


Rpq  =  |^a2  +  [qdy  -  (y  -  p  y ) 2  ]  ^ J 


2  ^  ,  .  -.2 
z  +  (pdx  -  x) 


(A- 14) 


(A-15) 


The  equation  for  the  array  vector  potential,  (A-13),  can  be  transformed  to  a  more 
convenient  form  using  Fourier  transform  methods.  The  required  Fourier  transform 
pair  is 


F(w)  = 


e 


(A-16) 


-•y[a2  +  (w-w1)2]1’ 

[a2  +  (w-w^2]^ 

jwxt  _ 

f(t)  =  — — j —  Hq2^  (a/-y2-t2)  where  w^  is  a  constant  (A-17) 

The  function  Hq2\x)  is  the  0^  order  Hankel  function  of  the  second  kind  and  f(t) 
and  f (w)  are  related  by  the  usual  Fourier  transform  equations. 


F  (w) 


f(t)dt 


(A-18) 


f(t) 


F  (w)  dw 


(A-19) 


Using  (A-16)  through  (A-19) ,  the  array  row  vector  potential  becomes 


_°°  .  - j kqd  s 

e  y  y 


F(qdy) 


(A-20) 


or,  more  conveniently. 


jqV 

e  F(qwQ) 


where 


(A-21) 


t  =  -ks 

y 

The  infinite  series  in  (A-21)  now  can  be  transformed  to  a  faster  converging  series 
using  the  Poisson  sum  rule  (shifted  version)17  which  states 


ejqv 

q*-co 


F(qwQ) 


f(t  +  nT) 


(A-22) 


A-4 


for  any  Fourier  transform  pair  F(w)  and  f(t).  The  parameter  T  is  given  by 


T  ,  2-n  =  2 it 

wn  d 
0  y 

Using  (A-22),  the  array  row  vector  potential  in  (A-21)  is 


(A-23) 


dA  =  T 
P 


f(t  +  nT) 


(A-24) 


The  function  f(t  +  nT)  is  given  by 


jw-(t  +  nT) 


*(')-£ 


Hq2)  (a/^2  -  (t  +  nT)2) 


(A-25) 


where  the  constant  w^  is  picked  to  be 


w]_  =  y  -  pAy 

The  array  vector  potential  then  satisifies 


(A-26) 


dA  =  l 


^oodl 

dl 


~jkpdxSx  - 
e  X  X  dA 


(A-26) 


The  row  vector  potential  is 


dA 


jy[ks  + 

y  °y 


,  ,  2im 

-  jpAy  — 

y  u(2) 


H^'  (akr') 


(A-27) 


where 


:’  =  ti^]2  -  d-s  +  ~)2)h 


LjkJ 


y  d 


Y  = 


2u 


Interchange  of  summation  operations  and  change  of  sign  on  the  summation  index  in 
(A-26)  and  (A-27)  yields 


dA  =  l 


uI00d* 


4jd 


-jykfs  +  !*•) 

y  dA 


(A-28) 


y  n=-°° 


A  -5 


where 


dA  = 
n 


-jkpdx[sx  -  r^] 

J  X  y  H(2)  (akr’) 

0 


(A-29) 


where 


''  ’  [  'A>2  -  Isy  +  Sr’*  ] 


It  now  is  necessary  to  transform  (A-29)  using  the  Poisson  sum  rule.  To  begin  the 

(2) 

process,  the  Hankel  function,  Hq  (x) ,  is  written  in  terms  of  a  modified  Bessel  func¬ 
tion  of  the  second  kind,  K^(y). 

Vu> =  ij  ho2)<-^>  (a_3°) 

where  U  is  a  variable.  Using  (A- 30),  the  Hankel  function  in  (A-29)  becomes 


H<2)(akr')  =  Kq  |Wz2  +  (pdx  -  x)2J 


(A-31) 


where  the  variable  U  is  chosen  to  be 

°-/2  +  k2[’y+r]2 


The  row  vector  potential  in  (A-29)  then  becomes 


¥  K0  [Vy/z2  +  (pdx  -  X)2  ] 


(A-32) 


To  transform  (A-32)  using  the  Poisson  sum  rule,  the  following  Fourier  transform  pair 
is  required. 


g(t) 

G(w) 


^-K0[a  [B2  +  (t-C)2]*3  ] 
2  ?  h 

-B(ANw  )  -jwC 
z2  e _ 
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(A-33) 

(A- 34) 


where  A,  B,  and  C  are  constants.  The  row  vector  potential  in  (A-32)  takes  the  form 


dA  = 
n 


-JPTnfi 

e  8(P^0) 


t  =  PTr 


T0  =  dx 


“  *  k<»,  -  ¥£> 

X  y 

and  the  constants  A,  B,  and  ?  are  chosen  to  be 


A  =  U  =  A2  +  k2 


The  Poisson  sum  rule  for  this  transformation  takes  the  form 


"jpTon 


g(pTQ)  =  ~  Y'  G(mW  +  R) 


0  m=-<» 


W  = 

To 

Upon  transforming  (A-35)  using  (A-36)  the  row  vector  potential  becomes 


dA  = 
n  yd 


-jkx  sx 

9  L. 


mX  nAyX 


d  d  d 
x  x  y 


“V* 


where 


(A-35) 


(A-36) 


(A-37) 


The  total  array  vector  potential  then  becomes 


dA  =  t 


vlOO*1 
2yd  d 

x  y  n=-°°  m=-°° 


z  z 


-yr • g , 


(A-38) 


for  r • z  < 0  where 


84 


=  ^  [s  v  +  m  T —  [s  +  ^p-jy  ±  g  z 

’L'  \  VyJ  v  Ly  dyj  z 


X  x  yj 

The  array  vector  potential  in  (A-38)  thus  is  seen  to  be  representable  as  a  doubly 
infinite  sum  of  attenuated  plane  waves.  The  fields  are  obtainable  from  the  vector 
potential  in  (A-38)  by  use  of  the  relations 


dS  =  i  V  x  dA 

w 


dE  -  ~  V  x  dH 
jwe 


The  final  results  are  expressed  in  terms  of  the  E  field. 


dE  = 


nIoodt 

2d  d 
x  y 


OO  OO 

z  z 

n=-oo  m=-» 


-jkr-g. 


(A-39) 


for  r-z?0  where 


hi  =  7  (i  X  g+)  x  g+ 
*z 
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APPENDIX  B 


LENS  CELL  SIZE 

The  size  of  a  cell  Is  limited  mainly  by  module-to-module  phase  settings  that 
differ  significantly  from  those  of  linearly  progressive  phasing.  The  first  approxi¬ 
mation  analysis,  applied  separately  to  each  cell,  assumes  linear  progressive  phasing 
between  modules  in  a  cell  (Sections  3.2.1  and  3.2.6).  A  nonlinear  progressive  phas¬ 
ing  across  a  flat  lens  is  necessary  if  a  target  return  is  to  be  focused  on  the  radar 
receiver;  however,  this  deviation  from  linear  will  be  minimal  across  a  sufficiently 
small  cell.  A  typical  SBR  lens  cell  size  is  determined  here.  A  circular  wavefront 
segment  is  indicated  in  Figure  B-l. 


Figure  B-l.  Linear  Modeling  of  a  Circular  Arc 
The  maximum  deviation  from  linear  is 


=  Tr2  +  rn  21 %  _  r 

L  l2J  J 


From  (B-l),  the  segment  length  for  a  specified  error,  ip,  is 
d  =  2  [(r+<j/)2  -  r2]h 


(B-l) 


(B-2) 


This  length  represents  the  height  of  a  cell  of  containing  modules,  with  linear  pro¬ 
gressive  phasing  that  approximates  the  actual  phasing  to  within  an  accuracy  of 
iJj  wavelengths. 

The  cell  height  required  to  maintain  a  specified  maximum  phase  error  will 
vary  depending  upon  the  distance  of  the  cell  from  the  lens  center.  The  relation¬ 
ship  of  cell  size  relative  to  distance  from  the  lens  center  is  show  in  Figure  B-2. 


B-l 


_  Y  +  d/2 
c 


l-Y  -  d/2 
c 

—  d0/2 


-d0/2 


Figure  B-2.  Cell  Height  Relative  to  Distance  from  Lens  Center 

The  radius  of  curvature  of  the  wavefront  at  the  lens  increases  as  the  distance  from 
the  lens  center  increases.  The  cell  size,  therefore,  also  increases  for  the  same 
accuracy;  for  example,  for  a  iJj  =  0.1A  deviation,  and  feed  to  lens  separation  of 
r  =  444A,  the  length,  d,  at  the  center  of  the  lens  is  d  =  18. 8A.  If  the  lens  radius 
is  approximately  146 A,  then  the  wavefront  curvature  radius  at  the  edge  of  the  lens 
is  r  =  468A  and  d  =  20A.  In  general,  the  cell  height,  d,  is  calculated  as  follows 
with  reference  to  Figure  B-2. 


9  =  tan 

c 


o  ►  "I  Yc 

9  =  tan  — 

r 

c 


where 


d0  =  cell  height  at  lens  center 


r  =  radius  from  receiver  to  cell  center 
c 

Y  =  distance  from  lens  center  to  cell  center 
c 


Y  +  d/2  =  r  tan  (6  +  0  ) 
c  c  c 


Y  -  d/2  =  r  tan  (9  -  0  ) 
c  c  c 


Therefore, 


d  =  r^  [tan  (0  +  0 )  -  tan  (9  -  9C)J 


For  the  above  example,  the  cell  height  at  the  edge  of  the  lens  (6c  =  1.29°,  0  *  18°) 
is  d  =  22. 2A,  while  at  the  lens  center  d„  ■  18. 8A. 


A  simple  approximation  to  the  number  of  cells  needed  to  model  a  lens  to  within 
an  error  of  ip  is  obtained  by  setting  all  cell  heights  d  =  d0.  The  number  of  cells 
needed  to  model  the  lens  then  is  calculated  by  dividing  the  lens  into  concentric 
rings  and  subdividing  the  rings  into  cells.  The  width  of  a  cell  is  set  equal 
to  the  height  of  the  cell.  The  number  of  cells  per  ring  increases  for  rings 
farther  from  the  lens  center.  In  the  above  example,  since  the  lens  radius  is 
148 X,  the  number  of  rings  needed  to  cover  the  lens  is  eight.  The  number  of 
cells  per  ring  is  shown  below  (0  ring  is  the  inside  ring). 

Ring  Number  Number  of  cells 
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The  total  number  of  cells  needed  to  cover  the  lens  is  204.  A  lens  portion 
subdivided  into  cells,  is  shown  in  Figure  4-3. 


MISSION 

of 

Rome  Air  Development  Center 

RA VC  plant,  and  exec ate*  research,  development,  tut  and 
detected  acquisition  programs  In  support  of  Command,  Contact 
Communications  and  Intelligence  (C3I)  activities .  Technical 
and  engineering  support  within  areas  of  technical  competence 
Is  provided  to  ESP  Program  Offices  (P04l  and  other  ESO 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  contact,  sur¬ 
veillance  of  ground  and  aerospace  objects,  intelligence  data 
collection  and  handling,  information  system  technology, 
ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 


